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ABSTRACT 


The Naval Postgraduate School is currently in possession 
of software designed to perform a thermal analysis of 
electronic components. This software package incorporates a 
model builder which contains two programs whose primary 
function is to generate a thermal model. In its present 
configuration, the model builder requires an inordinate amount 
of time for data input and model verification. This thesis 
describes the development of a model builder designed 
specifically to reduce the time required to model the 


Substrate,epoxy and carrier layers of a microcircuit assembly. 
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I. INTRODUCTION 


The Naval Postgraduate School 1s currently in possession 
of software designed to perform a thermal analysis of 
electronic components. This software package incorporates a 
model builder which contains two programs whose primary 
function is to generate a thermal model, or input data file, 
to be read by the thermal analyzer program. The first program 
is considered to be a general model builder which is used in 
all model development stages as well as to modify an existing 
model. The second alternative was developed to generate a 
thermal model of a specific microcircuit geometry. 

The development of an accurate thermal model of an 
electrical component requires that the structure be subdivided 
into a large number of small but finite subvolumes. Each 
subvolume is assumed to be isothermal with the centroids, also 
called nodes, considered to be representative of the entire 
subvolume. The most difficult problem encountered in the 


development of a thermal model is the generation of m-node 
equations in vA-unknown temperatures where the nodes are 


connected by thermal conductances. As the desired accuracy of 
the thermal model increases, the number of required node 


equations becomes extremely large. Therefore, it is imperative 


that the design engineer have access to a model builder that 
will produce the thermal model in a reasonable period of time. 

In its present configuration, the thermal analysis 
software contains a model builder that generates the required 
node equations automatically. There is no question that the 
existing model builder programs have replaced the extremely 
laborious and time consuming process of generating the node 
equations by hand. However, they still require an inordinate 
amount of time for data input and model verification. 

This thesis describes the development of a model builder 
designed specifically to reduce the time required to model the 
substrate, epoxy, and carrier layers of a é microciretm: 
assembly. A typical microcircuit package configuration is 
shown in Figure 1. Figure 2a provides a horizontal interior 
illustration while Figure 2b displays the specific geometry to 
be modeled. All three layers may contain an equal number of 
nodes over their width. However, the carrier layer may contain 
a mounting ear on the front and rear surfaces. Additional 


cHaracteryStics to be disclssea in what follows are. 


1) The capability of working t1neEnglish on Sivunsieee 
2) The choice of four aspect ratios. 


3) The provision for up to 740 nodes depending on the 
existence of mounting surfaces (ears). 


4) The ability to input heat dissipation using several 
methods. 


5) The provision for six ambient temperatures. 


6) The provision for rapid, menu-driven data input. 


7) The automatic calculation of conductance values based on 
user input. 





Biguee 1. TO-5 configuration... (Courtesy of Honeywell, inc.) 
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Li THE REASON FOR THERMAL ANALYSIS 


Over the past several decades a trend of increasing 
sophistication and complexity has enveloped the electronics 
industry. This continuing advance in technology has greatly 
increased the reliability, capability, performance, and 
avallability of electronic systems. The escalating demand for 
further advances in all areas of electronics has presented 
engineers with an abundance of complex problems. 

One major area of concern is the continued development of 
advanced methods in the thermal control of multilayered 
structures. It is the responsibility of designers to ensure 
that electronic components operate efficiently and effectively 
throughout the specified thermal limits. Therefore, it is 
extremely important that design engineers have the capability 
to accurately and rapidly predict the temperature distribution 
on multilayered structures prior to prototype production. The 
overriding reasons for performing a precise thermal analysis 
are to increase component reliability, ensure proper material 
selection, ensure bias stabilization, and reduce or eliminate 


the possibility of catastrophic thermal failure.{Ref. 1] 


A RELIABILITY 

There is a predictable relationship between the operating 
temperature of electronic components and reliability [Ref. 2]. 
The materials used in the fabrication of components have 
temperature limitations. Should these temperature boundaries 
be exceeded, the physical and chemical properties of the 
material are altered and the device fails. Figure 3 displays 
the intimate relationship between failure rate and component 
operating temperature for some selected devices. Furthermore, 
it is an established fact that the reliability of an 
electronic component is inversely proportional to the junction 
or component temperature and is also directly linked to 
failure rates [Ref. 1}. 

Consider Figure 4 which illustrates the "bathtub" 
mortality curve with the failure rate of a particular 
component plotted against component age during operation 
within thermal limits. The high failure rate in the interval 
prior to t,, also known asgthe Dumgoin (period ais sonsideeen 
to be the result of poor quality control during the 
fabrication process.[Ref. 1] 

The area of highest concern is the interval» between twand 
t,. This period is considered to be the useful lite) ss ies 
with proper guality control, testing and burn-in procedures, 


t, 1S equal to zero. Failures that occur in this interval are 


b 


due to a variety of causes and are unpredictable.[Ref. 1] 
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Figure 3. Failure rate vs. Temperature for 
selecte devices [Ref. 2]. 
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Figure 4. "Bathtub" mortality curve showing variation of 
failure rate with component age [Ref. 1]. 


Failure also occurs gradually due to sustained operations 
within specified temperature limits. If all aspects of 
faorication have been performed correctly only a small 
Mesetion of components will have failed prior to reaching t_, 
also known as the wear-out period. Failures that occur in this 
period are due to the slow and never ending deterioration of 


materials. 


B. MATERIAL SELECTION 

The fabrication of electronic components results in the 
joining of several different materials. Consider Figure 2 
which depicts a typical semiconductor structure. When power is 
applied to the components, heat is dissipated to the substrate 
and subsequently to the carrier. The mechanical properties of 
these materials are all affected differently by changes in 
temperature. Opposing mechanical and chemical reactions due to 
environmental conditions and contaminants may result in 
component performance degradation or a reduction in useful 
life. Table 1 lists temperature related factors that may 
affect component performance. 

The primary objective in the selection of materials for 
the fabrication of an electronic assembly is to achieve the 
desired level of correlation within the finished product. As 
packaging densities increase thermal, mechanical, electrical, 
and chemical coupling becomes very strong. This high level of 


coupling can be both an advantage and a disadvantage. For 


example, a high level of correlation is desirable during 
fabrication to ensure an uniform product. However, in use, 
strong coupling is generally more desirable for moderate 
temperature deviations and weaker coupling is more desirable 
for large temperature deviations. In the case of large 
temperature deviations strong coupling may result in the 
catastrophic failure of many connected components while weak 
coupling may limit the number of failed components. Therefore, 
the strength of coupling between materials must be based on 
the type of failure most likely to occur and an accurate 
thermal analysis must supply this information to assist in 


proper material selection [Ref. 3]. 


TABLE 1. TEMPERATURE FACTORS [Ref. 4] 


Effect on Equipment Accelerating factors 


——_ a RR : 
































Increasing Loss of strength, reduced Lubricants, rubber parts, 
Temperature stiffness, reduced resonant plastics, corrosion, fatigue, 
frequency, softening, distortion, load intensity, and time 






aging, and creep duration 












Lubricants, rubber parts, 
plastics, and time duration 






Reducing Temperature Increased viscosity, increased 
stiffness, increased resonant 
frequency, brittleness, and 
reduced impact resistance 



















Temperature cycling, temperature 
range, unequal expansion 
coefficients, stress 
concentrations 

and lack of strain relief 


Thermal Expansion and Change 1n size and shape, 
Contraction buckling, cracking, distortion, 
and loosening 
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Ce BIAS STABILIZATION 

The first step in the design and implementation of a 
semiconductor device is to establish a stable and predictable 
electrical operating point. This procedure, known as bias 
Stabilization, attempts to determine a stable operating point 
lia t is virtually independent Of external component 
parameters. However, aS external parameters change, the 
operating point is directly affected. Therefore, a good bias 
design ensures that components will always operate within a 
certain range of their nominal value.[Ref. 5] 

Consider Figure 5 which displays a transistor connected in 
the common-emitter configuration. Suppose that a proposed 
operation requires a specific collector to emitter voltage 
eee the Circuit consists of a battery or some other source 
Bates Provides a bias voltage V.., the collector resistor R,, 


and the transistor. By Kirchoff's voltage law 


= Vogt Ripert, Viera= 20 a 


which results in a collector to emitter voltage of 


(2) 


Should the collector current be allowed to increase in excess 


cepcolerable inits, V. must decrease because Vv. and KR, are 


CE 


uae 





Figure 5. Schematic of transistor connected in common-emitter 
configuration (tet. 1 jx 


Ie 


fixed values. Therefore, it should be noted that, if a high 
Bnecelon temperature Causes an increase in I), Vi- can no 
longer be maintained at the level desired to perform the 
desired operation.{Ref. 1] 
1. Operating in the Forward Bias Region 
As an example of electronic component temperature 
dependence, consider a diode operating in the forward bias 


region. In the forward region the i-v relationship is closely 


approximated by 


pee se 


Meelis equation I. 1S a constant for a given diode at a given 
m=mperacure. The current If, is usually called the saturation 
current. However, another name for it is the scale current, 
mabe arises from the fact that i. 1s directly proportional to 
the cross-sectional area of the diode. Furthermore, it can be 


seen in Table 2 that, I. 1s a very strong function of 


S 
temperature. [{Ref. 5] 
ihe temperature relationship between i. and the 


mengvard current 1 1S derived from the voltage V_. This 


constant, called the thermal voltage, is given by 


q 


13 


where 


k= Boltzman's Comertatiems | seumle:.- aoa 
T = the absolute temperature, K 
q = the charge on theselectronm, |se0o- 10 —« 


Table 2 illustrates this relationship and emphasizes the need 
to accurately analyze a proposed assembly prior to 


fabrication. 


TABLE 2. TEMPERATURE DEPENDENCY OF I, ON 1 


FOR SELECTED MATERIALS 
GERMANIUM SILICON 


5 0...0+ 7% 62.293, 1 BA 
Bla . {ie 







D. CATASTROPHIC THERMAL FAILURE 

Another of the primary goals of techniques in advanced 
thermal control is to provide a thermal environment for a 
diversity of components that are in increasingly close 
proximity to each other. Figure 6 illustrates the increasing 
level of packaging densities. With increasing complexity comes 
an increased level of connections between dissimilar material 
and aegreater possibility for exceeding temperature 


limitations. Therefore, it is necessary that designers have 
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Pergure 6. State of the art in circuit complexity [Ref. 1]. 


ales. 


some knowledge of the possible range of component and 
environmental variations in order to prevent catastrophic 
thermal failure.[Ref. 1] 

Catastrophic thermal failure is defined as an immediate, 
thermally induced, total loss of electronic function in a 
specified component. This type of failure is the result of a 
component melting due to excessive temperature, a thermal 
fracture of the substrate or carrier, or a separation of leads 
and the external network. It 1S generally considered to be 
dependent on the local temperature field, operating history, 
and operating modes of the component. AS previously stated, a 
variety of problems arise when components are subjected to 
temperatures in excess of their rated limits. Furthermore, it 
1s extremely difficult to determine the precise temperature at 
which catastrophic failure may occur. The incorporation of an 
accurate thermal analysis, in combination with test and 
operating experience, may be used to generate a catastrophe 
free upper operating limit. These maximum allowable operating 
temperatures are used to generate the master thermal control 


configuration for the system. [Ref.1] 
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a Sls Oe HEAT TRANSFER 


Heat transfer is defined as all energy flows that arise as 
a result of temperature differences [{Ref. 6]. Because 
electronic components are not one hundred percent efficient, 
they produce heat as well as the desired output. In the case 
of semiconductor devices, heat develops in parts having low 
thermal efficiencies, such as the die. One of the major 
objectives of packaging is to develop an effective system for 
the removal of heat from these parts [Ref. 4]. It is 
imperative that design engineers understand all modes of heat 
transfer in order to incorporate an efficient method of heat 
removal into component designs. The modes of heat transfer are 


conduction, convection, and radiation. 


A. CONDUCTION 

Conduction is the transfer by molecular motion of heat 
between one part of a body to another part of the same body or 
between one body and another in physical contact [{Ref. 1]. 
Joseph Fourier, a French physicist, proposed that the rate of 
heat flow through a material by conduction is proportional to 
the area of the material normal to the heat flow path and to 


the temperature gradient along the heat flow path. 


le 


This proportionality is represented mathematically by 


(5) 
ar 

x =A 
q ax 


where the minus sign allows for a positive heat flow in the 
presence of a negative temperature gradient. The introduction 
of a proportionality constant, known as thermal conductivity, 
results in the following rate equation which describes this 


mechanism [Ref. 1]: 


a (6) 


where 


k thermal conductivity of the material, W/m-°C 


A = area of the heat flow path, m¢ 
aT/dx = change in temperature per unit length, °C/m 


q = rate of heat flow, W 


1. General Equation of Heat Conduction 
The first step in the analytical solution of a heat 
conduction problem for a given structure is to choose an 
orthogonal coordinate system such that the surfaces coincide 
with the boundary surfaces of the structure [Ref. 7]. in tHe 


case of the model builder developed in this thesis,the 
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rectangular coordinate system will be employed. The general 


equation of heat conduction is given as 


oe (7) 


Ox 


m(Kge) * =| =| a o(k-ge 


Ox Oy. OZ 


Then, assuming k, C, and p are independent of temperature, 


direction, and time, the resulting equation is 


2 Z 2 
tei Oe 8 GS 1 Or (8) 
Ox? oy? Oz? K a ot 
where 
T = tempéecactize.  C 
x, y, and z = cartesian coordinates, m 
t = time, sec 
k = thermal conductivity, W/m-°C 
q = internal heat generation, W/m 
oe = thermaledttrusivity «© k7 pe. m/sec 


There are several variations of the general equation 
of conduction. The first, known as the Fourier equation, 
provides a solution for a system that contains no heat 


sources: 











er, OT PT. 1 aT oe 


a2 


The second variation, known as the Poisson equation, 
supplies a solution for a system in which the temperature does 


not vary with time: 











Es 10 
ge (10) 


The third and final variation of the general equation 
of conduction provides a solution for a system void of heat 
sources and operating in steady state. The resulting equation, 
known as the Laplace equation, is given as 


CP. Cee Ges 


VE GEG (11) 
(op 4a oy? Ome 











2. Simple Plane Slab 
Consider Figure 7 which illustrates a simple plane 


slab with face temperatures T “and T,. Using only fone 


1 


dimension, equation (11) is reduced to 





By integrating twice and applying boundary conditions the 


temperature distribution across the slab is seen to be 


Tape en - T,) (12) 
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Figure 7. Conduction heat transfer through a simple 
plane slab. 


Zo 


Insertion of equation (12) into equation (6) produces 


solution for the heat flow across the slab: 


= ir, - T,) (13) 








Ohm's Law indicates a direct analogy between heat 
flow, equation (13), and the flow of electrical current 
through a resistor, V = RI. This electrothermal analogy is 
extremely useful in the solution of one dimensional, steady 
state problems without energy generation and will be developed 
further in the next section.[Ref. 8] 

3. Electrothermal Analog 

As previously stated, there is a direct analogy 

between heat flow across a simple plane slab, equation (13), 


and electrical current governed by Ohm's law: 


Seas 
he = (14) 


In this case, the analogous quantities are 


Current I + Heat Flow q 
Potential V + Temperature Difference AT 


Resistance R + Thermal Resistance R 
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It is easily seen that for the heat flow in a simple 
plane slab described by, equation (13), the thermal resistance 


Ts 


R = eae (15) 


I Es 
g kA 
The electrothermal analog for conduction across a simple plane 


slab is shown in Figure 8.[Ref. 1] 


B. CONVECTION 

Convection is defined as the process by which thermal 
energy is transferred to or from a solid by a fluid flowing 
past it. Should the fluid flow be the result of a temperature 
difference the phenomena is called natural or free convection. 
On the other hand, when a pump or fan causes the mass movement 
the process is called forced convection. [Ref. 1] 

Recall that at the interface between a solid and a fluid 
that heat is transferred by conduction and must obey Fourier's 
law, equation (6). Due to the difficulty encountered in 
accurately measuring the temperature gradient, Newton 
suggested that the surface heat transfer rate be related to 
the product of surface area and the temperature difference 
between the surface and the fluid. The results of this 


proposition lead to Newton's law of cooling: 


q = hACT, a Tf) oe) 


Zo 


Figure 8. 


Conduction Electrothermal Analog. 
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5 |e 


where i is a proportionality factor that has become known as 
the surface heat transfer coefficient.[(Ref. 1] 

From a comparison of Newton's law of cooling, equation 
(16), and Fourier's law, equation (6), one can derive that the 
surface heat transfer coefficient can be related to the 
thermal conductivity, the wall temperature gradient of the 
fluid, and the surface fluid temperature difference: 


fat 
ee ei = 45 (17) 


AAT Nae 


Therefore, any correlation between heat transfer coefficients 


must reflect the dependence of A on the thermal conductivity 


of the fluid and on the ratio of the wall temperature gradient 
to the temperature difference.[{Ref. 1] 
1. Electrothermal Analog 
The addition of heat transfer by convection to both 
surfaces of the simple plane slab of Figure 7 results in the 
configuration shown in Figure 9. In the convective case 


thermal resistance is represented by 


=m 
hA 
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Figure 9. Convective heat transfer on a simple plane slab. 


Therefore, the total thermal resistance is 


+ 








Alt All 
oe Eee 19 
h, . ‘| ae 





which is represented by the electrothermal analog shown in 
Figure 10. A simple consideration of circuit theory then shows 


iia tc: 


(P= ft) 
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C. RADIATION 

Heat transfer by radiation is the means by which thermal 
energy can be transmitted through a space without an 
intervening medium while obeying the laws of electromagnetics. 
Thermal radiation, while traveling at the speed of light, may 
be absorbed, reflected, or transmitted upon contact with a 
surface. An ideal black body absorbs all incident radiation 
and reflects and transmits none of it. The concept of the 
black body is useful because laws governing its radiation are 
Simple and many real bodies may be treated approximately as 
black bodies [Ref. 1]. 

Materials used in the fabrication of electronic components 


are classified as gray. Gray bodies are diffusely reflecting 


Cag 





R= 
ie 
1. a 
me E 
KA 
q 
R= i 
h.A 
2 
T, 
Figure Ox Gandue eiliem and convection electrothermal 
equivalent. 


28 


Opaque surfaces [Ref. 4]. These surfaces reflect equal amounts 
of energy over the thermal radiation spectrum (wavelengths of 
Boouc 0.1 pm to about 100 um) in all directions. The heat 
transfer efficiency of this mode depends on the configuration, 
the orientation, and the temperatures of the surfaces in the 
electronic assembly. 
1. Transformation of the General Radiation Equation 

The use of the thermal radiation equation in 
analytical studies is made difficult by its dependence on the 
fourth power relationship between the temperatures. Due to the 
nonlinear characteristics and the complexity of the 
calculations a computer program is the desired method to solve 
problems that have a significant transfer by radiation. The 


general equation for radiation interchange is: 


G = OF ,F_A (T, - Te) (21) 


where 


o = Stefan-Boltzman constant, 5.669:10% W/m?-> K’ 


F, = shape factor that accounts for the 
arrangement of the of the radiating source and 
absorbing receiver 


EP eeecticeivitly factor thar accounts for the 
ability of the source and receiver to emit or absorb 
radiation 


= 
ll 


temperature of the source, °K 


=. 
I 


temperature of the receiver, °K 


A = surface area, m* 


Zo 


It 1S important to note that the absolute temperature scale 
must be used when considering radiation. 

One way to handle computations involving heat transfer 
by radiation is to transform the general radiation equation, 
equation (21), into a configuration compatible with Fourier's 
law. In this case, linearization of the general radiation 
equation is the method used to produce the desired result. 
This is achieved by factoring the difference in the fourth 


power of the temperatures as follows: 


(To > te) = (fete) (i. re) 
= (re (eee em eG cies) (22) 


Inserting this into equatiens (2b eresules in 


G= OF Ate ene CT. = Fe) (23) 


i 


A radiative heat transfer coefficient may therefore be defined 


as 


Bis SoG FRAT. ee, Cae) (24a) 


OF 


h. = oF Pate Toe Tey (24b) 
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tien, substitutang fh into equation (23), radiation heat 
transfer may be treated exactly as convection at the boundary. 
A thermal resistance for radiation heat transfer can 


now be proposed: 


il 


rw: 


(25) 


Figure 11 provides an illustration of the electrothermal 
equivalent with the addition of radiation resistance in 


parallel with convective resistance. 














Ty 
R= — R= = 
hy A ha 
setae TH 
| i 
R= 
kA 
Ty» 
1 1 
R= R= 
hr A h,A 
T 


Figure 11. Radiation, convection, and conduction 
electrothermal equivalent. 
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2. General Problems of Heat Transfer by Radiation 

In the calculation of heat transfer by radiation, it 
is usually necessary to approximate real~-body behavior by the 
gray-body idealization. The assumption that emissivity is 
always equal to absorptivity is frequently reguired if the 
problem is to be solved. Even a Simple situation becomes quite 
complex 1f real body behavior 1s considered and the resulting 
convenience of replacing absorptivity with emissivity is lost. 
A second difficulty in considering real-body behavior is the 
lack of sufficient data. To properly account for real-body 
behavior, extensive tabulations would be required.[Ref. 6] 

The usual problem of heat transfer by radiation is 
further complicated not only because the surfaces are nonblack 
but also because the configuration of the areas involved is 
not simple and reflecting surfaces may be present to augment 
the direct exchange [Ref. 6]. There are a number of methods 
available to assist in the solution of heat transfer by 
radiation, however, their development 1s beyond the scope of 


this -thesa.. 


BZ 


ee FINITE DIFFERENCES 


The need to utilize a computer to determine the 
temperature distribution within an electronic component or 
system has made finite difference methods very desirable. This 
numerical method is useful lia problems involving 
nonlinearities, complex geometries, complicated boundary 
conditions, or a system of coupled partial differential 
equations. The purpose of this section is to provide the 
reader with some basic concepts involved in finite difference 
methods for solving differential equations. Furthermore, it 


demonstrates the methodology used to formulate nr#-node 
equations in m-unknown temperatures and instills confidence 


that this numerical method is capable of generating an 


accurate thermal model. 


A. FUNDAMENTAL CONCEPTS 

Consider equation (7), the general equation for heat 
transfer by conduction. In order to produce a numerical 
solution to a conduction heat transfer problem it is necessary 
to reconfigure the partial differential equation into a form 
that allows differentiation to be performed by numerical 
methods. Therefore, lee is essential that accurate 
approximations of the first and second derivatives be 


obtained. 
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1. First Derivative Approximation 
The derivative of a function at any point can be 
expressed by a finite difference approximation by 
incorporating a Taylor series expansion about that point. 
Consider Figure 12 where T(x) is a function that can be 
expanded by a Taylor series. The Taylor series expansions of 


the functions T(x+Ax) and T(x-Ax) about a point x are: 


3 
T(x+Ax) = T(x) + axt’(x) + 18% pcx) « AO mE +... (26a) 


Tix-Ax) = FX) - axr’(x) + 182 p(x) - S49 r(x) ae (26b) 


In order to determine the first derivative, equations (26a) 


and (26b) are solved for T’ (x). 


T(x) = Dhe+ Ax) - TOO _ (AX) py yy ~ (AXP gy -) (27) 
AX 2 

1 (Xe) aed oe Oe) ee eae (28) 
AT Z 6 


Subtracting equation (26b) from (26a) and solving for Tag 


produces 


Ti(x) = T(x + Ax) - T(x - Ax) _ (AX)? p(y) Se (29) 
2AX 6 
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First and second derivative approximations. 


Figure 12. 
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From equations (27), (28), and (29) the fiteu 


derivative approximations of T’(x) are determined to be 


T(x) = = = Ti ag Ti) forward difference (30) 
T(x) = = = Thx) = AM) backward difference (31) 
Tix) = a = Ev ax Pa axt central difference (32) 


and for all it is observed that 


2. Second Derivative Approximation 
To obtain the second derivative of the function T(x), 
Figure 12, consider the Taylor series expansions of the 


functions T(x+2Ax) and T(x-2Ax) about a point x. 


T(x+2Ax) = T(x) + 2AxT’(x) + 2(Ax)?T"(x) + = (Ax) 37""(x) ee (33a) 


TUx-28x) = TCs) = PAT Aare = (Ax)? 7 (x) Stal (33b) 
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Inserting the corresponding first derivative approximation 


into equations (33a) and (33b) and solving for T” (x) results 








in 
T(x) = Jat SES og PEO Sa eA. NS a 6 ee et ee ee (34) 
ox (Ax)? 
p(x) = LE = Ble > 2hx) + Tx) ~ 2M ~ Ax) backward difference (35) 
ax CA 


The central difference 1s obtained by eliminating T’ (x) 


between equations (26a) and (26b). 





TH! (x) = ash = EU CSL SMEs GS Me. 9 central! difference (36) 
ax¢ (Ax) 2 
and for all it is observed that 
Z 
TH! (x) = omy = Ca) 
ax (Aga 


B. NODE ANALYSIS 

The first step in the physical formulation of a solution 
by node analysis 1S, as previously stated, to divide the 
region into a finite number of subvolumes as shown in 
Figure 13. The centroid of each subvolume is called a node and 


is considered to be representative of the entire subvolume. 
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Bigure 12.) SErucEuTe Givided aneo 


At this point, continuity must be applied. In the steady 
state, the algebraic sum of all the heat leaving any node must 
equal zero. As an example, consider an interior node on the 
upper layer; select node number 14. Figure 14 provides an 
expanded picture around node-14. In order to maintain clarity 
the external environment will be considered as a single node, 
number 81. The resulting energy balance or node equation can 


be written as 


Oi ion 2 Gin een oe ia aie 8 aes. * Siar * Gz; = O (37) 


Here, each numerical subscript represents the heat flow from 
node-14 to the indicated node. Furthermore, q, allows for 
energy input by an external source such as a dissipating 
component. 

The energy balance simply supplies the pertinent energy 
terms. It is important to note that differences in thicknesses 
and thermal conductivities must be taken into consideration. 
The rate equations may be written in terms of thermal 


resistances and temperature differences: 


k, AxAZ 

Q14,13 > ae Ty 3) (38a) 
k,AxAz 

W14,15 ~ : Ay - Tie = 5) (38b) 
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where 








M4,10 = aa By) (38c) 

Qi4a,e1 = een, = are) (38e) 

cee aa, BE, Ca) oe, 
K, k, 


kK andes thermal Condquctivities of the top layer 
and middle layer, respectively, W/m-°C 


Ax and Ay = path lengths in their respective 
direction and are node dimensions for surface 
area calculations, m 


AZ Zon teAZ ee sene thicknesses Tor the upper and 
middle layers, respectively, m 


T = temperature of the node indicated by the 


subscript, ce 


Substituting these rate equations into the energy balance 


equation, equation (37), results in an equation containing 


five unknown temperatures. 
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For simplicity let Ax—Ay-AZ ja2)=1 and kk keen 


K(T,, — Tyo) + K(Ty, - Ty3) + K(Ty, -— Tys) + K(T,, — Tye) 


+ 2 = eae ae) ope (39) 


firtiher, reductions fresh hes 1a 


-Ti9 ~ Tig > Tis = eee eee eee a Seen (40) 


An equation of this type can be written for every node in the 
region whose temperature is unknown. This results ina set of 


nm-equations that relate the n-unknown temperatures. There are 


a variety of methods available to solve for the unknown 
temperatures. In large structures the most desirable method 


would be a thermal analysis program. 
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V. THE MODEL BUILDER 


As previously stated, the model builder currently being 
used by the thermal analysis software package requires an 
inordinate amount of time for data entry. This thesis presents 
a specific model builder, TMDL, that will become part of the 
software package. TMDL provides a user friendly, rapid method 
to assist in the development of a thermal model with a 
specific geometry as in Figure 2. Specifically, it generates 
a properly formatted output data file for use by the thermal 
analyzer enroute to preparing an accurate listing of the 


temperature distribution of the structure. 


A. THE THERMAL ANALYZER INPUT DATA FILE 

The model builder generates a data file from physical 
characteristics of the structure provided by the user. This 
data file, also called the thermal analyzer input data file, 
must be in a format that is completely acceptable to the 
thermal analyzer. Furthermore, it will consist of five lines 
and up to seven data sets. This section describes each line 
and data set and their relationship with TMDL. Figure 15 
displays a partial data file. 

Line one is the title line. It may be left blank or may 
contain up to 79 alphanumeric characters. The user selected 


title appears at the top of the data file. 
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EXAMPLE OUTPUT DATA FILE 


740 6 0 0 0 0 0 0 2 
0 0 0 

750 50 6 2 4 6 0 0 0 

0 

-0500000 .6666700 12 .8000000 78.00000 


76.500000000 76.900000000 75.900000000 77.100000000 77.400000000 75.6000 


7 JooL 21 7521 111 7511 2411 
134.100 36.850 87.709 43.855 2.047 1.456 
6 11 31 7521 121 7511 2421 
36.850 36.850 87.709 43.855 2.047 1.456 
7 21 41 7521 sl Toit 2431 
36.850 36.850 87.709 43.855 2.047 1.456 
6 31 51 1521 141 7511 2441 
36.850 36.850 87.709 43.855 2.047 1.456 
6 41 61 7521 151 #511 2451 
36.850 36.850 87.709 43.855 2.047 1.456 
6 a1 71 7521 161 g511 2461 
36.850 36.850 87.709 43.855 2.047 1.456 
7 61 81 7521 171 7511 2471 
36.850 36.850 87.709 43.855 2.047 1.456 
6 gl 91 7521 181 7511 2481 
36.850 36.850 87.709 43.855 2.047 1.456 
6 81 101 7521 191 Howl 2491 
36.850 36.850 87.709 43.855 2.047 1.456 
q 91 7541 V28 201 F51t 2501 
36.850 73.700 87.709 43.855 2.047 1.456 
6 7551 121 11 2.1 Tore 2511 
73.700 36.850 43.855 43.855 2.047 1.456 
6 111 131 21 221 751% 2521 
36.850 36.850 43.855 43.855 2.047 1.456 
7 121 141 31 231 7511 2531 
36.850 36.850 43.855 43.855 2.047 1.456 
6 131 151 41 241 7511 2541 
36.850 36.850 43.855 43.855 2.047 1.456 
6 141 161 51 251 7511 2552 
36.850 36.850 43.855 43.855 2.047 1.456 
6 151 a7 61 261 J511 2561 
36.850 36.850 43.855 43.855 2.047 1.456 
6 161 181 Lek 271 Tol 2571 
36.850 36.850 43.855 43.855 2.047 1.456 
6 171 191 81 281 7511 2581 
36.850 36.850 43.855 43.855 2.047 1.456 
7 181 201 91 291 7511 259) 
36.850 36.850 43.855 43.855 2.047 1.456 
6 191 7541 101 gO 7511 2601 
36.850 73.700 43.855 43.855 2.047 1.456 
6 7551 221 111 311 Toil 2611 
73.700 36.850 43.855 43.855 2.047 1.456 
6 211 231 121 321 Lore 2621 
36.850 36.850 43.855 43.855 2.047 1.456 
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Figure 15. Output data file. 


Line two is the problem data line. It has nine entries of 
which two are under user control, number of nodes’ under 
consideration and unit type. One entry, the number of constant 
temperatures is preset at six for this specific model. The 
remaining entries have applications to models associated with 
heaters, unique exponents, secondary heat input, temperature 
coefficients and curves, and nodes controlling fast heat. 
These entries are not applicable to this model and are preset 
to zero. 

Line three places a zero at three points and 1s beyond the 
user's control. Therefore, no further discussion is required. 

Line four is the problem capability line. This line 
defines the maximum values for the entries in line two. The 
first entry is 750, the number of nodes for which the analysis 
is dimensioned. This 1S significant because the first constant 
temperature will be assigned the node number 751. The second 
entry 1s 50 which is the maximum number of constant 
temperatures in accordance with the analyzer dimension 
statement. Therefore, it 1S possible to have 50 constant 
temperatures allocated from node 751 to 800. The third entry 
is preset to six for an application concerning heaters and is 
not applicable to this model. The balance of the entries in 
line four represent a listing of data sets that are required 
for the particular analysis at hand. TMDL uses three data sets 


miatc will be discussed in what follows. 
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Line five contains five items that concern the level of 
accuracy that the thermal analyzer will achieve. These entries 
are preset. Item one provides the desired level of accuracy 
between iterations. Because the thermal analyzer solves 
iteratively, there must be an error criterion at which point 
calculations cease. A number that is too small will cause the 
computer to run excessively and a number that is too large 
will not give the desired accuracy. Therefore, a trade off is 
necessary and such a tradeoff has indicated that 0.05 is 
Satisfactory. Item two is a damping factor that is used 
between iterations to prevent temperature oscillations between 
iterations. The third entry is the maximum number of 
iterations. This prevents excessive computer time in the event 
of faulty input data. Item four is a convergence factor that 
adjusts the damping to close to the critical value. This means 
that once the computer determines convergence is occurring, 
slow convergence is increased to reduce computer time. The 
fifth and final entry in line five is the initial temperature 
at which the iterative process begins. It is input by the 
user. 

Input data set one contains temperature dependent 
coefficients and is not used in this model. 

Input data set two contains up to 50 constant temperature 
inputs. This model has six constant ambient temperatures input 


by the user. 
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Input data set three involves heaters for fast warm up and 
is not used in this model. 
Input data set four contains all pertinent information 


concerning the w-node equations. Each node requires two lines 


of data. The odd numbered lines are used for specifying the 
nodes that interact with the node in question and the modes by 
which this interaction takes place. Consider Figure 14, line 
one of data set four is as follows; 

6 Ggeiayl a Ges yo dt ad 7 Sa: . <2 
This line of data is for node number one of aspect ratio 
selection one. The first entry provides the number of 
connections to that node. Entry two says node 755, an ambient 
temperature, 1S connected to the node in question and the one 
indicates that the connection is by conduction. The rest of 
the entries are read similarly with numbers between 7511 and 
7561 indicating a conductive connection with the external 
environment. External heat input is defined by entry 9991. 

The even lines are used for specifying the inter node 
conductance values of the corresponding entry in the previous 
line. 

Input data set five is used only if there are unique 
exponents. Therefore, no further discussion is required. 

Input data set six contains the initial temperature 


guesses corresponding to the number of nodes’ receiving 
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secondary heat and is not used in this model because they have 
been set by the last entry in line five. 
Input data set seven indicates the number of temperature 


dependent heat input curves and is not used in this model. 


DB FEATURES 

TMDL is presented to the user in discrete sections. Upon 
entry into the first section, the user is offered an optional 
overview. This option should be accepted on at least the first 
run. At the completion of the overview the user is prompted 
for a data file name and title. The data file name should be 
changed for each successive run because TMDL does not over 
write existing files. 

The second section provides prompts for the physical 
characteristics of the structure. Initially, the user is given 


a list of four aspect ratios from which to choose; 


1.) 10 by 24 nodes provides a 1:2.4 ratio with 720 nodes. 
2.) 15 by 15 nodes provides a 1:1 ratio with 675 nodes. 
3.) 12 by 20 nodes provides a 1:1.6 ratio with 720 nodes. 


4.) 8 by 30 nodes provides a 1:3.75 ratio with 720 nodes. 


It should be noted that up to 30 additional nodes may be added 
depending on the existence of the mounting ear. Figure 16 
displays the specific geometry with definitions. 

After selecting the desired aspect ratio, the user must 


specify the use of either SI or English units. The input of 
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Figure 16. 
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data must remain consistent with this selection. At this point 
the input of structure characteristics are requested. These 
entries consist of the length, width, thickness, and thermal 
conductivity of each layer. There is ample opportunity for 
correction or alteration throughout TMDL. 

Section three requests the input of initial and ambient 
temperatures. TMDL has s1x ambient temperatures. 

The final section allows for the input of external heat 
sources. TMDL provides for heat input into the upper substrate 
surface only. There are four methods of external heat input 
from which to choose. The first alternative allows for a total 
rate of heat applied to the upper surface. An entry for this 
alternative would be divided by the nodes and distributed 
appropriately. Option two provides for the entry of average 
heat per unit area. The third alternative gives the user the 
freedom to enter heat in specifically selected nodes. 

Option four offers no heat input. 

At the conclusion of these entries an output data file 

will be created in the proper format and placed in a file 


named by the user. 
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APPENDIX A 


TMDL LISTING 


SLARGE 

c 

CewILlE: MODEL BUILDER 

Gee 0HOR: LY PATRIC ROESCH 
CeDATE: 15 MARCI 1991 

C COMPILER: MICROSOFT VERSION 4.01 
@eeeINKER: MICROSOFT VERSION 3.55 
C 

© DEFINE REAL VARIABLES 


eee OW DEE She i Obri KE Ce Cw DELCKCITUPRI,.LWRI,LILRT 

cies THEAT HPN ArIEAI NEAT DELXDELYSYLR.SXFBSAX,SYY,SZE,SZT 
Pee ees Pie Yo or 7G lee Y PRE. CY YE. CXFBE.CXXE,CZEAR.CYLR,CYY,CXX 
ime2b,ACC. DAMP, CONFAC EYLR 


Qe ke 


DEFINE INTEGER VARIABLES 
INTEGER NPL,NWIDE,NDEEP,NUM,NUMA,CH,H, TOTNOD,NN,NC,N IB JE,IC,1D, 
+COUNT,CONTEMP,USEL,ZER,NMAX,TMAX,HTRS,D1,D2,D3,D5,D4,D6,D7,MAXIT 


2) 


DEFINE ALL ONE CHARACTER VARIABLES 
CHARACTER*1 ANS,ANSN,ANSA,ANSSL,ANSSW,ANSST,ANSSK,ANSE,ANSEL, 
+ ANSEW,ANSET,ANSEK,ANSC,ANSCL,ANSCW,ANSCT,ANSCK,ANST,ANSTLANSTU, 
+ ANSTLR,ANSTL,ANSTR,ANS TF,ANSTB,ANSH,ATH,ANSHA,AIIN,SLLECT,SSELH 
tANSWER,DELT 


CHARACTER VARIABLES OF MORL THAN ONE POSITION 
CHARACTER DATAF*70,UK*11,UT*1,UH* 13, UAH*6,NAME*6,UL*2 


DEFINE MATRICES 
REAL HEAT(30,15),COEF(740,9) 


INTEGER IH(240),JH(240),NCON(740,9) 


C 
C 
C 
C 
@ 
c 
C 
é 
C PROVIDE THE USER WITH A PROGRAM OVERVIEW 
'e 
1006 CALL CLS 

WRITE(*,1001) 
1001. FORMAT(///,’ oh 

' THIS PROGRAM WAS WRITTEN TO INTEGRATE WITH EXISTING | ’/, 

THERMAL ANALYSIS SOFTWARE AND TO REDUCE THE AMOUNT _ ’/, 
OF TIME REQUIRED FOR DATA ENTRY.’/, 
> WOULD YOU LIKE AN OVERVIEW OF THE PROGRAM PRIOR TO _ ’/, 
’ BEGINNING? ENTER Y FOR YES OR N FOR NO. ’,2X,\) 

READ(*,1002) OVR 
1002. FORMAT(A1) 


+++++ 
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G 
IF (OVR.EQ.’Y’.OR.OVR.EQ.’y’) THEN 
GADINGIEs 
WRITE(*,1003) 
1003 FORMAT(///,’ PAYEE ER Vib Wwoett tthe We 


+’ THIS PROGRAM PERFORMS A NODAL ANALYSIS OF A THREE __ ’/, 
+’ LAYERED CONDUCTIVE STRUCTURE (SUBSTRATE,EPOXY,AND  ’y, 
+’ CARRIER). THE OUTPUT CONSISTS OF UP TO 740 COEFFI- */, 
+’ CIENTS THAT CONTRIBUTE TO THE DETERMINATION OF  ‘y, 
+’ THE TEMPERATURE DISTRIBUTION OF A MICROCIRCUIT */, 
+° WHEN FED INTO THE THERMAL ANALYZER. ‘/, 
+’ THE FOLLOWING IS AN ORDERED OUTLINE OF THE MAJOR _’y, 
+’ SECTIONS OF THIS PROGRAM AND WHAT ENTRIES ARE __ */, 
+’ REQUIRED OF THE USER. _’/, 
+’ NOTE THAT ALL ENTRIES WILL BE IN UPPER CASE LETTERS. ’/, 
Yo SADA TAUE IEE es 
a 1.) ENTER THE TITLE OF THE OUTPUT DATA FILE. ‘/, 
ie 2.) ENTER THE DATA FILE NAME. THIS PROGRAM WILL NOT’, 
if ERASE OR WRITE OVER EXISTING FILES. °// 
+’ PLEASE PRESS ENTER WHEN READY TO CONTINUE. °,2X,\) 
READ(*,1002)ANSWER 
C 
CATING 
WRITE(*, 1004) 
1004 ==FORMAT(///, u 


+B.) STRUCTURE PHYSICAL CHARACTERISTICS 

Me 1.) SELECT UNIT TYPE (SI OR ENGLISH) _ ’/, 

o 2.) SELECT FROM FOUR ALTERNATIVES THE DESIRED NODAL’/, 

ae ASPECT RATIO. ’/, 

rs 3.) ENTER EACH LAYER CHARACTERISTICS. SUBSTRATE IS ’/, 

- THE TOP LAYER, FOLLOWED BY EPOXY AND THE CARRIER’, 

+ LAYERS, RESPECTIVELY. THE CARRIER LAYER MAY _ ’/, 

- CONTAIN AN EAR ON THE FRONT AND BACK SURFACES ’/, 

ss IF SPECIFICATIONS REQUIRE IT.’/, 

+’ PLEASE PRESS ENTER WHEN READY TO CONTINUE. °.2X,\) 
READ(*,1002)ANSWER 


CALL CLS 
WRITE(*, 1005) 
1005 = FORMAT(///, 
C.) INITIAL AND AMBIENT TEMPERATURES _ ’/, 
1.) ENTER INITIAL CHIP AND CHASIS TEMPERATURE. ‘/, 
ASSUME STEADY STATE FOR INITIAL TEMPERATURE. ’/, 
2.) ENTER AMBIENT TEMPERATURES FROM ALL SIDES OF THE’/, 
STRUCTURE.’,/, 
D.) HEAT INPUT ’/, 
HEAT INJECTION OCCURS ONLY ON THE UPPER SURFACE OF ’/, 
THE SUBSTRATE LAYER. THIS PROGRAM SUPPLIES THE USER ’/, 
FOR ALTERNATIVE METHODS FOR ENTERING HEAT; ’/, 
1.) TOTAL HEAT OVER SUFACE _ °/, 
2.) AVERAGE HEAT PER UNIT AREA ’/, 
3.) INPUT HEAT NODE BY NODE ’°/, 
4.) NO HEAT INPUT ’/, 
THIS COMPLETES THE PROGRAM OVERVIEW. PLEASE PRESS ENTER ’/, 
TO CONTINUE. °,2X,\) 


t+etttttt+ttetttetett 
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READ(*,1002)ANSWER 
C 
ELSEIF (OVR.EQ.’N’.OR.OVR.EQ.’n’) THEN 
GOTO 1007 
ELSE 
GOTO 1006 
1007 ENDIF 


c 
C 
es 
CeeNITIALIZE MATRICES 
Cc 

DATA HEAT /450*0.0/ 
DATA IIT /240*0/ 
DATA JH /240*0/ 


VARIABLE, CONSTANT, AND STRING DEFINITION 
PHYSICAL CHARACTERISTICS 
SL,EL,CL - SUBSTRATE,EPOXY, AND CARRIER LENGTHS 
SW.EW,CW - SUBSTRATE,EPOXY, AND CARRIER WIDTHS 
KS,KE,KC - SUBSRATE,EPOXY,AND CARRIFR THERMAL CONDUCTIVITIES 
DELS,DELE,DELC - SUBSTRATE,EPOXY AND CARRIER THICKNESS 
DELX - SL/NDEEP 
DELY - SW/NWIDE 
NPL - NUMBER OF NODES PER A LAYER 
NWIDE - NUMBER OF NODES WIDE 
NDEEP - NUMBER OF NODES DEEP 
UL - UNITS OF LENGTH (SI OR ENGLISH) 
UK - UNITS OF THERMAL CONDUCTIVITY (SI OR ENGLISH) 
EAR - DEPTH OF THE EAR (SL-CL)/2 


INITIAL AND AMBIENT TEMPERATURES 
IT - INITIAL CHASIS AND CHIP TEMPERATURE 
LT - LEFT SIDE AMBIENT TEMPERATURE 
RT - RIGHT SIDE AMBIENT TEMPERATURE 
FT - FRONT AMBIENT TEMPERATURE 
BT - BACK AMBIENT TEMPERATURE 
UPRT - UPPER AMBIENT TEMPERATURE 
LWRT - LOWER AMBIENT TEMPERATURE 
UT - UNITS OF TEMPERATURE (CENTIGRADE OR FARENHIET) 


meal INPUT 

(MHieEAT - TOTAL INJECTED HEAT 

taPN = TOTAL IEAT FER NODE 

AHEAT - AVERAGE HEAT OVER A GIVEN SURFACE 

NHEAT-=HEAT PER NODE INJECTED NODE BY NODE 

UH AND UAH - UNITS OF HEAT (SI OR ENGLISH) 

NUM,NUMA,CH,H - DUMMY VARIABLES USED TO CREATE THE VECTORS 
Niece 1 WHICH ARE USED TO RELATE NODE 
NUMBER TO MATRIX POSITION 

TOTNOD,NN,NC.N - VARIABLES USED TO ALLOW HEAT INPUT NODALLY 

IH. JH - VECTORS USED TO CORRELATE NODE NUMBER WITH 

MATRIX POSITION 
HEAT - MATRIX USED TO CONTAIN HEAT INPUTS 


CHO? OY OO 1) OO Or OO OO) OC Cee e@ a ere Oe ee Cree 2G a) a 
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PiQiovere elo tplere ore.e alevele Groreve erene) eve one etoler are: her enovee! ene ele@Lehee @ cel elelele eee. 


COEFFICIENT DEFINITIONS. YIMPLIES WIDTH & XX IMPEIES DER i 
ZIMPMIES IEG 
SIMPEMES SUBSTRATE RE IMREIESSELONY 
C IMPLIES CARRIER 


SYLR = PROVIDES COEFFICIENT TOR CE OR RIGHTEDGE NODES 
TO THE EXTERNAL NODE. SYLR IMPUEIES SUBSTKATE. 
EYLR = EPO xy 
CYER - CARRIER NO EAR 
CYDPRE -CARKIEK w/ EAK 


SXFB = PROVIDES COEFFICIENT FOR FRONT AND BACK EDGE N@DES 
TO THE EXPER NAL NODE: 
EXPE - EPO, 
CXFB - CARRIER/ NO EAR 
CXFBE - CARRIER W/ EAR 


SYY - PROVIDES INTERNAL COEFRIGIENT IN DHE Y DIRECTION 
EYY =EPO <1 
CYY = CARRIER W/NO EA 
CY YE = CARRIBRe Ww) bAK 


SXX= PROVIDE INTERNAL COEFFICIENT IN THE 'X DIREGHGN 
EXX - EPOXY 
CXX - CARRIER W/NO EAR 
CXXE - CARRIER W/ EAR 


SZ1T - COEFFICIENT FOR SUBS IRATE TOCUPECR EXTERNAL NODE 
SZE + COEPBICIENT FOR SUBSTRATE WOsErOx, 

EZC - COEFFICIENT FOR EPO ey TO GAkRiIER 

CZEAR = COEFFICIENT FOR EAR TOwlo. ENGExXTEKNAL NODE 
CZB - COEFFICIENT FOR CARRIER TO TON EREXTEKNAL NODE 


CHARACTER 
ALL ONE CHARACTER STRINGS EXCEPT Ul ARE SIMPLE 
YES, NO, OR SELECTION ANSWERS AND DO NOT REOUIE 
EXPLANATION. 


DATA FILE 

DATAP = TITLE FOR DATA FILE, DINE ONE 

NAME - DATA FILE NAME 

IB,IE,ID,IC - COUNTERS FOR FILLING DATA FILE 

COUNT - TOTAL NUMBER OF NODES 

ACC - MINIMUM ERROR CRITERIA AT WHICH CALCULATIONS CEASE 

DAMP - DAMPING FACTOR FOR PREVENTION OF TEMPERATURE 
OSCILLATIONS 

CONFAC - CONVERGENCE FACTOR USED TO INCREASE DAMPING 

CONTEMP - TMDL HAS 6 CONSTANT AMBIENT TEMPERATURES 

USEL - DELESANALYZER WHAT Vy Pr ORIN Ts evn IN ws 

ZER - DUMMY VARIABLE 

NMAX - MAXIMUM NODES PROGRAM IS CAPABLE OF EVALUATING 

TMAX - MAXIMUM NUMBER OF CONSTANT TEMPERATURES PROGRAM 
CAN ACCEPT 

HTRS - NOT APPLICABLE TO THIS PROGRAM - NO EXPLANATION 
NECCESSARY 
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& D1-7 - DATA SETS IN USE 
ec MAXIT - MAXIMUM ITERATIONS COMPUTER WILL PERFORM TO 
eo ACHIEVE SOLUTION 
Cc 
C BEGIN PROGRAM 
e 
Gx CLS 
WRITE(*,400) 
406 FORMAT (///,’ aR AE ok ook oe AB oe Ee os fs of oft os sR sR oS EB oR os oS ok oR Ook Bok oh ok ok 
res” | 
+ PRs ee ee a aa Se ee ee Ae ES 
ae AK sd 
+! -* PREASE-SEERGIGUPVER- CASE PERTERS ee. 
+° _ PRIOR TO BEGINNING. aa 
“e, ok ee. 
a ee ee ee Oe EE Se eS 
+’ ne a aan ge eee ee ie See 
Cc 
Se (*,201) 
201 format(//,’ ye 
+’ PRIOR TO ENTERING DATA INTO THIS PROGRAM ENSURE THAT’, 
+’ YOU HAVE A DRAWING OF YOUR DESIGN AND ALL PERTINENT ’,/, 
+’ DATA. ie 
+’ PRESS ENTER WitEN READY TO CONTINUEZ2ZX\) 
READ(*,304)ANSWER 
C 
GeeeieOW USER TO NAME THE DATA FILE 
GALE CLS 
Cc 
WRITE(*,701) 


701 FORMAT(///,, THIS PROGRAM CREATES AN OUTPUT DATA FILE FOR ENTRY 
+ INTO THE’), 
+’ EXISTING THERMAL ANALYZER. FURTHERMORE, THIS PROGRAM DOES’ /, 
+* NOT ERASE OR WRITE OVER THE EXISTING DATA FILE. THEREFORE’, 
+’ THE USER WILL NAME THE DATA FILE FOR EACH RUN OF THIS ’/, 
+’ PROGRAM. THE FILE NAME IS LIMITED TO SIX CHARACTERS. © ‘//, 
+” PLEASE ENTER THE DESIRED DATA FILF NAME. ’,2X,\) 
READ 702, NAME 

702 FORMAT(A6) 

C 

€ 

c 

C ALLOW USER TO PROVIDE TITLE LINE FOR DATA FILE 
GALL CLS 
WRITE(*,202) 

202 FORMAT(////,’ ENTER THE DESIRED TITLE 1O BE PLACED ON LINE 

+’ /. NUMBER-ONE‘ OF THE OUTPUT DATA FILE. /, 

+’ go) 

READ 1,DATAF 

FORMAT(A70) 


GAEL CLS 


COR Gis 


SURELY A LISMNG- OF ACCEPTABLE ASPECT RATIOS 
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C 
WRITE(*,203) 
203 FORMAT(/////,, THE FOLLOWING ASPECT RATIOS ARE AVAILABLE.’ 
WRITE(*,204) 
204 FORMAT(/,” 1.) 10 BY 24 PROVIDES A 1:2.4 RATIO WITH 740 NODES. 
+’/” 2.) 15 BY 15 PROVIDES A 1:1 RATIO WITH 675 NODES.’/, 
+’ 3.) 12 BY 20 PROVIDES A 1:1.6 RATIO WITH 720 NODES.’,/ 
+, 4.) 8 BY 30 PROVIDES A 1:3.75 RATIO WITH 720 NODES.’,// 
+, PLEASE SELECT A NUMBER ONE THROUGH FOUR. °,2X,\) 
READ(*,3020)SELECT 
020 FORMAT(A1) 


(2 


Ge® 


USER CAN CONTINUE OR MAKE ANOTHER SELECTION 
CAMInGEs 
IF (SELECT.NE. 1.AND.SELECT.NL.’ 2. AND.SELECT.NE.'3.AND.SELECT. 
+NE.’4’) THEN 
GOTO 7 
ELSE 
10  WRITE(*,3) SELECT 
3 FORMAT(////,; YOU SELECTED NUMBER ’,A1,’ OF THE FOLLO 
+WING 4 ALTERNATIVES.’,/, 
e 1.) 10 BY 24 FOR A 1:2.4 RATIO’, 
ze 2.) 15 BY 15 FOR A 1:1 RATIO’, 
a 3.) 12 BY 20 FOR A 1:1.6 RATIO’/, 
- 4.) 8 BY 30 FOR A 1:3.75 RATIO‘//) 
e 
ENDIF 
@ 
5 WRITE(*,205) 
205 FORMAT(’ IS THIS THE DESIRED SELECTION. ENTER Y FOR YES AN 
+D N FOR NO. °,2X\\) 
READ(*, 8)ANS 
FORMAT(A1) 


er 


IF (ANS.EQ.’N’) THEN 
GOTO 7 

ELSEIF (ANS.EQ.’Y’) THEN 
GOTO 9 

ELSE 
Atle lcs 
GOTO 10 

ENDIF 


fo ee 


CALIZCKES 


Ee 

é 

C AFTER CONFIRMATION OF SELECTION FILL CONSTANTS WITH APPROPRIATE 

C VALUES 
IF (SELECT.EQ.’1’) THEN 


NWIDE = 10 
NDEEP = 24 
NPL = 240 
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ELSEIF (SELECT.:Q.'2’) THLN 


NWIDE = 15 
NREEP = 15 
Nem 225 


Por lP (SELECT IEO.3') TEEN 


NWIDE = 12 
NDEFP = 20 
NPL = 240 
C 
ELSEIF (SELECT.EQ.'4’) THEN 
NWIDE = 8 
NDEEFP = 30 
NPL = 240 
ENDIF 
C 
C MAKE DESIRED UNIT SELECTION 
Cc 
2 WRITE(*,206) 
206 FORMAT(//////, | THIS PROGRAM IS CAPABLE OF OPERATIONS IN EI 
+THER SI OR ‘/, 
+’ ENGLISH UNITS. AFTER THE SELECTION OF THE UNITS, ALL ’/, 
+) ENTRIES MUST BE COMPATIBLE. PLEASE MAKE YOUR SELECTION.’/ 
+/, |S FOR SI NOTATION’, 
+’ EFOR ENGLISH NOTATION ’,2X,\) 


READ(*,302)ANSN 
302. FORMAT(A1) 
EC 
C CHECK FOR CORRECT UNFIT SELECTION 
4 IF (ANSN.EQ.’S’) THEN 
WRITE(*,207) 
207  FORMAT(////,, |= YOU SELECTED SI NOTATION.’) 
ELSEIF (ANSN.EQ.’E’) THEN 
WRITE(*,208) 
208  FORMAT(////’ ©§ YOU HAVE SELECTED ENGLISH NOTATION.’) 
ELSE 
GALL CLS 
GOTO 2 
ENDIF 
C 
WRITE(*,209) 
209 FORMAT(’) — IS THIS THE DESIRED SELECTION? ENTER Y FOR YES AND 
+N FOR NO. °,2X,\) 


READ(*,303)ANSA 
303. FORMAT(A1) 
@ 
C 


C SHOULD I CONTINUE OR RESELECT UNITS 
IF (ANSA.EQ.’Y’) THEN 
GOTO 6 
ELSEIF (ANSA.EQ.'N’) TIIEN 
GA CLS 
GOTO 2 
ELSE 
eA Ges 
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WRITE(*,210) 
210  FORMAT(//////, THIS PROGRAM IS CAPABLE OF OPERATIONS IN 
+ EITHER SI OR’/, 
+’ ENGLISH UNITS. AFTER TIIE SELECTION OF UNITS. ALL’/, 
+’ ENTRIES MUST BE COMPATIBLE. PLEASE MAKE YOUR SELECTION.’/ 
+/,  § FOR SI NOTATION.’ 
+’ E FOR ENGLISH NOTATION.’ 


G 
GOrO4 
6 ENDIF 
C 
CALL-CrES 
Cc 


(CARO OIG ie ioioi Gi i GR Sok aca acacia: ak cae SG i rae a eae CAC ee HOR RC Oke Re ae aC ok fe oe coe oe cai ae aie fs 2 a a 


Ctthtsseennaeaees® ENTER CHIP CHARACTERISMES HOR EACHELANmeIe ttt tts aan 
(COO GI RG ao Gi IGG a icici ag aoe ie ac ac a ae cok Ae co oie at coke eo oka Ake oe ee eo ede eee oie ac 

C 

C 


€ 
WRITE(*,211) 


a | FORMAT (fre ee ee eee ee! a te eee ae 
> ’ 
ee He ae a he ae a a ea le 2 zie 
4 OR A kai ke EAE CTI RSTRATE CHARACTERISTICS * *%% "#44 8% 84 4% #% 
jee Ree 
ae DIOR IGIICICIOIOROROIOIC I ROIGIOIOI IOI IG IIR IGIGIO ISRO IOI OIOIORCROIOROIC ICICI TOI IG IOIOIOI CHC HCHO HCHO HCHO Hee 


se oie ae ak oe ake a ie? 
+ /) 


C PROVIDE CORKECT UNIT ABBREVIATIONS 
IPZANSN-EO: SS) bEIEN 
WRIPE(~ 212) 
as FORMAT( ALL ENTRIES ARE IN SI NOTATION.’/ 
UR= sen 
UK = ’Watts/em/C 
8 eae O 
BCSEIF CANSN bOSE THEN 
WIE C213) 
213 FORMATC ALL ENTRIES ARE IN ENGLISH NOTATION.’,/) 
ie =a ih 
UR =. Biushr/in/l’ 
OS oes 
ENDIF 


WRITE(*,220)UL 
220 FORMAT(’ ENTER SUBSTRATE LENGTH (’,A2,’): a, 
READ *,SL 


WRITE(*,216)UL 

216 FORMAT(/, ENTER SUBSTRATE WIDTH (’,A2,’): 2a 
READ *,SW 
WRITE(*,217)UL 

217. FORMAT(/, ENTER SUBSTRATE THICKNESS (’,A2,’): _’,2X,\) 
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READ * DELS 
c 
hy Rie (* 218)UK 
218 FORMAT(/,, ENTER SUBSTRATE TIIERMAL CONDUCTIVITY (’,A12,’): °,2X,\) 
READ *,KS 
2 
ieee ALIL CLS 
eS 
C MAKE ANY CHANGES OR CORRECTIONS TO SUBSRATE ENTRIES 
© 
WRITE(*,214) 
214 FORMAT (///,’ YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE SUB 
POT RALE.’/,) 
Pome r RINT *, 1.) LENGTH ’SL 
PRINT *,’ 2.) WIDE “SW 
BRIN *, a TNC RNESS SPELLS 
PRINT *, +) k “KS 
eC 
Rie (*,215) 
215 FORMAT(, DO YOU WISH TO MAKE ANY CHANGES? SELECT Y FOR YES 
+ AND N FOR NO. ’,2X,\) 
READ(*,304)ANSS 
304. FORMAT(A1) 
C 
ie (ANSS.EQ.’Y') THEN 
12 CALL: CES 
WY RITE(*,219) 
219 FORMAT(//// 
PRINT * THE CURRENT ENTRY FOR LENGTH IS SL,’ ’,UL 
WRITE(*,2190) 
Zee rORMAT(/, WOULD YOU LIKE TO CHANGE THE LENGTH? (Y OR N) 
oes, \) 
READ(*,304)ANSSL 
PRINT * 
IW ANSSUCEO,Y ) THEN 
WRITE 221) UL 
221 FORME sete Niel. SUBSTRATE LENGIH.(A2.): *.2X,\) 
READ *,SL 
BiSEIF (ANSSL.EOUN) TEEN 
GOTO 14 
ELSE 
GOTO 12 
ENDIF 


© 
14 CALL CLS 
WRITE(*,2191) 
2191 FORMAT(////) 
PRINT *, Ine CURRENT ENTRY FOR WIDTH IS Sw, ’,UL 
WRITE(*,222) 
Zoe FORMAT, SevOUurs YOU LIKE TO CHANGE THE WIDTH? (Y OR N) |, 
+2X,\) 
READ(*,2192) ANSSW 
2192 FORMAT(AI) 
IF (ANSSW.EQ.'Y’) THEN 
WRITE(*,223)UL 
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223 FORMAT(/, ENTER SUBSTRATE WIDTH (’,A2,’): ’.2X,\) 
READ *,SW 
ELSFIF (ANSSW.LQ.N’) THEN 
GOTO 13 
ELSE 
GOTO 14 
ENDIF 


13 GAIEE-GLS 
WRITE G2193) 
2193 FORMAT (///) 
PRINT *, THE CURRENT ENTRY FOR THICKNESS IS DELS UL 
WRITE(*,224) 
224 FORMAT(/,’ WOULD YOU LIKE TO CHANGE THE AHIGKNESS: 41 ORay 
+ ’,2X,\) 
READ(*,2194) ANSST 
2194 FORMAT(AL) 
IF(ANSS VEO! ¥*) THEN 
WRITE 22a et 
225 FORMAT(/, ENTER SUBSTRATE THI@KNESS (AZ. 
SEA) 
KEAD * DELS 
EESEIe {ANSS EON ) THEN 


GOTO 15 
EESE 
GOTO 13 
ENDIF 
Cc 
€ 


15 CALL GIES 
WRITE(*,2195) 
2195 FORMAT(///) 
PRINT *7 SHE CURRENT ENTRY FOR THERMAL CONDUCTIVITY IS4ks 
+, OK 
WRITE(*226) 
226 FORMAT(/, WOULD YOU LIKE TO CHANGE THE THERMAL CONDUCHIVIT 
+Y? (XOORIN) 2 2X\) 
READ(*,2196) ANSSK 
2196 FORMAT(A1) 
IF (ANSSK.EO.’Y’) THEN 
WRITE 227)UK 
may FORMAT(/, ENTER SUBSTRATE THERMAL CONDUCTIVITY (,A1 
+2,’): ’,2X,\) 
READ *,KS 
ELSEIF (ANSSK.EQ.’N’) THEN 
GOTO 16 
EESE 
GOTO 15 
ENDIF 
Cc 
C ALLOW ANOTHER REVIEW OF SUBSTRATE DATA ENTRIES 
C 
16 CALL CLS 
WRITEGE228) 
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228 FORMAT(///,’ YOU HAVE MADE TIFE FOLLOWING CORRECTIONS TO T 
Pile SUBSTRATE ENTRIES.) 


GOTO 20 
ce 
Cc 
Pieieite (ANSS le. N ) 1 MIGIN 
re 10) 17 
Piste 
GOO 31 
ENDIF 
c 
cS 


CREE EAEKEREEEKEKK EARS ARH AEREE DON Y CELA RACTERIS TICS *** # #4 * # # * # # # #4 4 # & & 
C 
@ 
Pow CALL CLS 

WRITE(*,229) 
229 Be ey ee se toe ee ee et ee ee eee eee 

ee EN ee ee? 

oF 4 


EE EM lak lai aia aa dialed) oA OP. CHARAC TERIS [ICs itt ttt rrr es tree 


KK ke 
ct 8 


= Ea malay | 


fae te PAAT SUBSTRATE AND EPOXY LENGTH AND WIDTH ARE EQUAL. 7 
+) 


WRITE(*,230)UL 

230 FORMAT(/, ENTER EPOXY THICKNESS (’,A2,’): ’,2X,\) 
READ *,DELE 
WRITE(*,231)UK 

231 FORMAT(/, ENTER EPOXY THERMAL CONDUCTIVITY (.A12,): ,2X,\) 
READ *,KE 


C 
Bi— SL 
EW = SW 
G 
eee IEW EPOXY ENTRIES 
C 


30 CALL CLS 
WRITE(*,232) 

232. FORMAT(///, YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE EPOXY 
+ LAYER.’/) 


32. PRINT *, 1.) LENGTH EL 
PRINT *,’ 2.) WIDTH EW 
PRINT *, 3.) THICKNESS ’ DELE 
PRINT *, 4.) k KE 

C 


WRITE(*,233) 

233. FORMAT(/, NOTE THAT CHANGING EPOXY LENGTH OR WIDTH ALSO CHANG 
+ES’/, SUBSTRATE LENGTH OR WIDTH.’,/, 
+’ DO YOU WISH TO MAKE ANY CHANGES TO THE EPOXY ENTRIES? (Y OR 
Ne °2X.\) 
READ(*,2197) ANSE 

2197 FORMAT(A1) 

C 


61 


C MAKE ANY CHANGES OR CORRECTIONS TO EPOXY ENTRIES 
C 
IF (ANSE.EQ.’Y’) THEN 
22 CAERLE CLS 
WRITE(*,2198) 
2198 FORMAT(///) 
PRINT? THE CURRENT ENIRY FOR VENGIBIS 2 cr 
WRITE(*,234) 
234 FORMAT(/, WOULD YOU LIKE [TO CHANGE THE LENG EH Yy-OKaN) 
+2) 
READ (*,2199) ANSEL 
2199 FORMAT(AI) 


IF (ANSEL.FQ.Y°) THEN 
WRITE(*,235)UL 
235 FORMAT(/, ENTER EPOXY LENGTH (’,A2.’): ’,2X,\) 
READ *,EL 
SL = EL 
ELSEIF (ANSEL.EQ.’N’) THEN 
GOTO 21 
ELSE 
GOTO 2? 
21 ENDIF 


24 GARIRGES 
WRITE(*,601) 
601 FORMAT(////) 
PRINT ™*, THE CURRENT ENTRY FORSMADIT ISs.EW, {UE 
WRITE(*,236) 
236 FORMAT(/,, WOULD YOU LIKE TO CHANGE THE WIDTH? (Y OR N) 
+’,2X,\) 


READ(*,602)ANSEW 
602. FORMAT(A1) 
€ 
IF (ANSEW.EQ.’Y’) THEN 
WRITE(*,237)UL 
237 FORMAT(/, ENTER EPOXY WIDTH (’,A2,): ‘.2X,\) 
READ *,EW 
SW = EW 
ELSEIF (ANSEW.EQ.'N’) THEN 
GOTO 23 
ELSE 
GOTO 24 
23 ENDIF 
e 
2 CALLCLS 
WRITE(*,603) 


603 FORMAT(////) 
PRINT *, THEsGURREND ENTRY EOR THICKNESS IS. .DELE, UL 
WRITE(*,238) 

238 FORMAT(/,’ WOULD YOU LIKE TO CHANGE THE THICKNESS? (Y OR N) 

+ ’,2X,\) 

READ(*,604)ANSET 

604 = FORMAT(A1) 

& 
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ie (ANSETEO, Y') THN 
WRITE(*,239)UL 

239 FORMAT(/. ENTER EPOXY PHICKNESS (Az. ): |.2X,\) 

READ? DELI 
PESEM (ANSETEO.N’?) TEEN 

GOEO 25 
BESE 

GOTO 26 
ENDIF 


ho —e) 
Ce 


CAEL CLS 
WRITE(*,605) 
605 FORMAT (///) 
PRINT *,) THE CURRENT ENTRY FOR THERMAL CONDUCTIVITY IS ’,KE 
+,” (UK 
WRITE(*,240) 
240 FORMAT (/,’ WOULD YOU LIKE TO CHANGE THE TIERMAL CONDUCTIVIT 
oe yY OR N) .2X,\) 
READ (*,606)ANSEK 
606 FORMAT(AL) 


c 
nT ANSEKO.Y ) THEN 
WRITE(*,241)UK 
241 FORM AG UNTER EPOXY lHERMAL-CONDUCIIVITY (,Alz2,) 
on 2X,\) 
Eo ele (ANSEK.EO.N') THEN 
COTO, 
ELSE 
GOTO 28 
ENDIF 
C 
Gee OVIDE ANOTHER REVIEW OF EPOXY ENTRIES 
C 
Pac ALL CLS 
RITE (* 242) 


242 FORMAT(//, YOU HAVE MADE THE FOLLOWING CORRECTIONS TO THE EP 
+OXY ENTRIES.’//) 


GOTO 32 
E 
ELSEIF (ANSE.EQ.’N’) THEN 
GOTO 29 
ELSE 
GOTO 30 
ENDIF 
@ 
E 
et PEASE EERE REAR EKER AR RTE TR CHARACTERISTICS © ttt oe ter te tee et 
c 
e 
PomeALL CLS 
WRITE(*,243) 


43 FORMAT(//// Ok AR I aa GG aOR IC a Ga a ak ae ak de a0 dt ake fio dct deb of at ake akc ote dade of de ake ode ote te ake ok 
s. 


ae a a oR as ae ook ok a ak oo Ak ak ak a? 
cf He 


LEEK NAAR AAA eee AR RIER CHARACTERISTICS (erent ste hee 
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aK ag aK ok oe AS KS 
a i, 


te a ord) 

E 

944 WRITE(*,244)UL 

244. FORMAT(’) ENTER CARRIER LENGTHE (,A2,): .2X,\) 
READ *.CL 


ae 
C CARRIER LENGITIT MUST BE GREATER THAN OR EQUAL TO SUBSTRATE LENGTH 
S 


IF (CL.LT.SL) THEN 
CALL CLS 
WRITE(*,243) 
PRINT *,, CARRIER LENGTH MUST BE GREATER THAN OR EQUAL TO SUBS 
+TRATE LENGTH. ° 
GOTO 944 
ENDIF 
ie 
C 
WRITE(*,246)UL 
246 FORMAT(/,, ENTER CARRIER THICKNESS (’,A2,"): ’,2X,\) 
READ *,DELC 
C 
WRITE(*,247)UK 
247. FORMAT(/,, EN’TER CARRIER THERMAL CONDUCTIVITY (',A12,’): ,2X,\) 


REARS nC 
Cc 
33 GCAEICkS 
CW = SW 


C 
CVASEVIEW CAKKIEI@IN TISIES 
o 


WRITE(*,248) 
248 FORMAT(///”, YOU HAVE MADE THE FOLLOWING ENTRIES FOR THE CARRI 
+ER.’/) 


36 PRINT *, 1.) LENGTH Hele 
PRINT *, 2.) WIDTH "Cw 
PRINT *, 3.) THICKNESS = ’\DELC 
PRINT *, 4.) k *Kie 
PRINT * 


PRINT *,). CHANGING WIDTH WILL ALSO CHANGE SUBSTRATE AND EPOXY W 
+IDTHS.’ 
C 
WRITE(*,249) 
249 FORMAT(/, DO YOU WISH TO MAKE ANY CHANGES? (Y OR N) ',2X,\) 
READ(*,607)ANSC 
607 FORMAT(A1) 
E 
C MAKE CORRECTIONS OR CHANGES TO CARRIER ENTRIES 
€ 
IF (ANSC.EQ."Y’) FHEN 
33. FG@ANEIECIts 
WRITE(*,608) 
608  FORMAT(///) 
933 PRINT *, TITIE CURRENT ENTRY FOR LENGTH IS ’,CL,’ ,UL 
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WRITE(*,250) 
250 FORMAT(/, WOULD YOU LIKE TO CHANGE LENGTH? (Y OR N) °,2X, 
+\) 
READ(*,609)ANSCL 
609  FORMAT(AI) 


C 
IF (ANSCL.EQ.’Y’) THEN 
WRITE(*,251)UL 
251 FORMAT(/’ ENTER CARRIER ILENGTH (’,A2,’):.2X,\) 
READ *,CL 
HF (CL.LI.SL) THEN 
CALL CLS 
WRITE(*,010) 
610 FORMAT(////) 
PRINT *,, CARRIER LENGTH MUST BE GREATER THAN OR EQ 
+UAL TO SUBSTRATE LENGTH.’ 
PRINT *,” THI: PRESENT ENTRY FOR SUBSTRATE LENGTH IS 
est Ul, 
PRINT * 
GOTO 933 
ENDIF 
eC 
ELSEIF (ANSCL.EQ.’N’) THEN 
IF (CL.LT.SL) THEN 
CALL CLS 
WRITE(*,611) 
611 FORMAT(////) 
PRINT *,’ CARRIER LENGTH MUST BE GREATER THAN OR EQ 
+UAL TO SUBSTRATE LENGTH. ’ 
PRINT *” THE PRESENT ENTRY FOR SUBSTRATE LENGTH IS 
ae UL 
PRINT * 
GOTO 933 
ENDIF 
GOTO 39 
ELSE 
GOTO 33 
ENDIF 
C 


39 =9CALL CLS 
WRITE(*,612) 

612. FORMAT(////) 
PRINT *, THE CURRENT ENTRY FOR WIDTH IS CW,’ UL 
PRINT *,, CHANGING THIS ENTRY WILL CHANGE SUBSTRATE AND ’ 
PRINT *,’ EPOXY WIDTHS. THEY ARE ALL EQUAL. ’ 


WRITE(*,252) 
252. FORMAT(/, WOULD YOU LIKE TO CHANGE WIDTH? (Y OR N) °,2X,\ 
i) 
READ(*,613)ANSCW 
613. FORMAT(A1) 
c 
IF (ANSCW.EQ.’Y’) THEN 
WRITE(*,253)UL 
253 FORMAT(/” ENTER CARRIER WIDTH (’,A2,’):’,2X,\) 
READ *,CW 
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SW = CW 
EW = CW 

ELSEIF (ANSCW.EQ.’N’) THEN 
GOTO 34 

ELSE 
GOTO 39 

ENDIF 


oe) 
roe ®, 


CARL Cis 
WRITE(*,614) 
614 FORMAT(////) 
FRINT =. « tHE CURRENT ENTRY FOR THICKNESS, DEEG Ur 
WRITE(C* 25%) 
254 FORMAT(?, WOULD YOU LIKE TO CHANGE THICKNESS? (Wek Nj) 
+o 
READ(*,615)ANSCT 
615 FORMAT(A]1) 
G 
LE.CaSCleO” Yo) TREN 
WRI 2a hUL 
259) FORMAT,  ENREREGAICR TIP EIGKNESS (\A2,): ‘.2X,\) 
READ =, BELG 
ELSEIF (ANSCILEO.N ) TaEN 
GOTO'35 
EBSie 
GOTO 34 
ENDIF 


55 GAIT CLS 
WRITE(* 616) 
616 FORMAT(////) 
PRINT *” THE CURRENT ENTRY FOR THERMAL CONDUCTIVITY IS#7ke6 
+ UK 
WRITE(*,256) 
256 FORMAT(/”, WOULD YOU LIKE TO CHANGE THE THERMAL CONDUCTIVIT 
tY2 (Y OR Nj ogee 


READ(*,617)ANSCK 
617 FORMAT(A1) 
'@ 
IF (ANSCK.EOY’) THEN 
WRITE(@ 25)UK 
a) FORMAT(/”) = ENTER CARRIER THERMAL CONDUCTIVITY (,A12, 
+ 2X) 
READ * KC 
EESEIE {ANSCK.EO. NN?) THEN 
GOTO 40 
ELSE 
GOTO'35 
ENDIF 
C 
C ALLOW FOR ANOTHER REVIEW OF CARRIER ENTRIES 
C 
40 CALL CLS 
WRITE(*,258) 


258 FORMAT(////” YOU HAVE MADE THE FOLLOWING CORRECTIONS TO THE C 
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+ARRIER ENTRIES.’//) 
GOTO 36 


ELSEIF (ANSC.EQ.’N’) THEN 
Goro 37 
ELSE 
GOTO 38 
37. ENDIF 
C 
C 
& 


(Ce a ao ae etek ek a Ack eR TEN PERATURE [NPUTS *# ##% #4 8 #8 8 4 8 8 ee a ke ea kok ik ee a 
(UH adc ccc aGioi diicidaididick  aiggiciioiidciicigdgiicd icici doicicicde ddd acd a aca acdsee at at ea 
C 

CALL CLS 

WRITE(*,260) 


260 CIS TR TOT a aaa i al a a al tal 
eee! | 
cf EE EE EEE See ANMIBTIEN | TEMPERATURE Nie Woe 
ae 
ma cedcns aston ee WE HHH EAE He A AE EAN NE GR 6 8 DE NE i 2 Se I eS 
oe eff) 

C 

C INITIAL TEMPERATURE 

Cc 
mRITE(*,.261)0UT 

201 FORMAT(  LNTER THE INITIAL CHIP TEMPERATURE (’,Al,’). °,2X,\) 
READ *,IT 

Ee 

C UPPER AMBIENT FEMPERATURE 

e 
WRITE(*,262)UT 

262 SFORMAT(/, ENTER THE UPPER SURFACE AMBIENT TEMPERATURE (’,Al,’ 
=) ,2X,\) 


READ * UPRT 
Cc 
C LOWER AMBIENT TEMPERATURE 
C 
eel PE(*,2603)UT 
Poo ORMAT(?/, ENTER THE LOWER SURFACE AMBIENT TEMPERATURE (,Al,’ 
me 2X.\) 
READ *,LWRT 
G 
C RIGHT SIDE AMBIENT TEMPERATURE 
C 
WRITE(*,264)UT 
264 FORMAT(/, ENTER THE RIGHT SURFACE AMBIENT TEMPERATURE (’,Al,’ 
cle «2 s.Nh 
READ * i. Ii 
C 
@eerl SIDE AMBIENT TEMPERATURE 
G 


WRITE(*,265)UT 
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205 FORMAT(/, ENTER THE LEFT SURFACE AMBIENT TEMPERATURE (Ar 
+. ',2X,\) 
READ it 

C 

C FRONT AMBIENT TEMPERATURE 

e 
WRITEC 266) U1, 

206 FORMAT(/° ENTER THE FRONT SURFACE AMBIENT TEMPERATURE (’,Al,’ 
eee) 
iA lan 

C 

C REAR SIDE AMBIENT TEMPERATURE 

C 
WRITE(*,267)UT 

267 FORMAT(/, ENTER THE REAR SURFACE AMBIENT TEMPERATURE (’,A1,’) 
Fee aoe 
READ BT 


REVIEW TEMPERATURE ENTRIES 


Ge) Ge 


CALL CLS 

WRITE(*,268) 

208 FORMAT(////,, YOU HAVE MADE THE FOLLOWING AMBIENT TEMPERATURE 
+ENTRIES.’/) 


ERENT 1.) [INDIALSCHIP TEMPERATUIS “DE UT 
PRINT *,” 2.) UPPER AMBIENT TEMPERATURE OIA, Ua 
PRINT *, 3.) LOWER AMBIENT TEMPERATURE iy Re 
PRINT *, 4.) RIGHT AMBIENT TEMPERATURE wRT, UT 
PRINT *, 5.) LEFT AMBIENT TEMPERATURE pd ee nn) B® 
PRINT *, 6.) FRONT AMBIENT TEMPERATURE le) 
PRINT *, 7.) REAR AMBIENT TEMPERATURE eal 

CG 

e 


WRITE(*,269) 
209 FORMAT(/,,. WOULD YOU LIKE TO MAKE ANY CORRECTIONS? (Y OR N) ’ 
+,2X,\) 
READ(*,618)ANST 
618 FORMAT(A1) 
C 
C MAKE CORRECTIONS OR CHANGES TO TEMPERATURE ENTRIES 
E 
IF (ANST.EQ.’Y’) THEN 
42 \ GANIC ES 
WRITE(*,268) 


PRESS *, 1.) INITIAL CHIP TEMPERATUKE mae: {Om 

PRINT =, 2.) UPPER AMBIENT TEMPERATURE VEL Ur 
PRINT *,” 3.) LOWER AMBIENT TEMPERATURE SLWRI, UT 
PRINT *,’ 4.) RIGHT AMBIENT TEMPERATURE Ae UT 
PRINT, 5.) ERT AMBIENT TEMPERATURE alle UT 
PRINT *,” 6.) FRONT AMBIENT TEMPERATURE ol AO HY 
PRINT *, 7.) REAR AMBIENT TEMPERATURE (Bly UT 


c 
WRITE(*,2682) 
2682. =FORMAT(///, WHICH TEMPERATURE WOULD YOU LIKE TO CHANGE?’ / 
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+.) SELECT ZERO THROUGH SEVEN. ,2X,\) 
READ(*,618)DELT 
e 
C CORRECT INITIAL TEMPERATURE 
IF (DELT.EQ.’1’) THEN 
@.0 CLs 
WRITE(*,270) 
270 FORMAT(///) 
PRINT *” INITIAL CHIP TEMPREATURE IS “IT, ‘\UT 
WRITE(*,271)UT 
ae FORMAT(/” ENTER THE NEW VALUE,(,A1,)_ ’,2X,\) 
READ *,IT 


C CORRECT UPPER AMBIENT TEMPERATURE 
ELSEIF (DELT.EQ.’2’) THEN 
GALL CLS 
WRITE (*,620) 
620 FORMAT (//) 
PRINT *, UPPER AMBIENT TEMPERATURE IS ’,UPRT, ’,UT 
WRITE(*,.273)UT 
273 FORMAT(/, ENTER THE NEW VALUE.(,A1,) ’,2X,\) 
READ *,UPRT 


C CORRECT LOWER AMBIENT TEMPERATURE 
ELSEIF (DELT.EQ.’3') THEN 
CALL CLS 
WRITE(*,623) 
623 FORMAT(/) 
PRINT *,, LOWER AMBIENT TEMPERATURE IS ’,_LWRT,’ UT 
WRITE(*,275)UT 
275 FORMAT(/, ENTER THE NEW VALUE.(’,A1,’) ’,2X,\) 
READ *,LWRT 
é 
C CORRECT RIGHT AMBIENT TEMPERATURE 
ELSEIF (DELT.EQ.’4’) THEN 
GALL CLS 
WRITE(*,626) 
626 FORMAT(///) 
PRINT * RIGHT AMBIENT TEMPERATURE IS ’\RT,’ UT 
WRITE(*,277)UT 
277 FORMAT(/, ENTER THE NEW VALUE.(,A1,’) \2X,\) 
READ *,RT 


C CORRECT LEFT AMBIENT TEMPERATURE 
ELSEIF (DELT.EQ.’S’) THEN 
CONICS 
WRITE(*,629) 
629 FORMAT(///) 
PRINT *,, LEFT AMBIENT TEMPERATURE IS ‘LT, ‘UT 
WRITE(*,279)UT 
279 FORMAT(/, ENTER THE NEW VALUE (’,A1,’) °,2X,\) 
READ *,LT 
C 
C CORRECT FRONT AMBIENT TEMPERATURE 
ELSEIF (DELT.EQ.’6’) THEN 
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CULES 
WRITE(*,632) 

FORMAT(///) 
PRINT * FRONT AMBIENT TEMPERATURE IS \PTY UT 
WRITE(*,281)UT 

FORMAT(/” ENTER THE NEW VALUE,(’,Al,’) ‘,2X,\) 
READ *,FT 


C CORRECT REAR AMBIENT TEMPERATURE 


ELSEIF (DELT.EQ.’7’) THEN 

CAliGles 
WRITE(*,635) 

FORMAT(///) 
PRINT *,,. REAR AMBIENT TEMPERATURE IS ’,BT,’ ’,UT 
WRITE(*,284) UT 

FORMAT(/, ENTER THE NEW VALUE,(’,A1,)_’,2X,\) 
READ *,BT 


C MAKE NO CHANGES 


ELSEIF (DELT.EQ.’0’) THEN 
CAIINGIS 
WRITE(*,6320) 

FORMAT(///,, YOU HAVE DECIDED TO MAKE NO CORRECTIONS!’/, 


+ iS Mls CORRECT 1 CRN) 2) 


o 


READ(*,618) ANSTI 

IF (ANSTLEQ.’Y’) THEN 
GOTO 51 

ELSEIF (ANSTLEQ.’N’) THEN 
GOTO 42 

ELSE 
GOTO 4201 

ENDIF 


ELSE 
GOTO 42 
ENDIF 


C ALLOW FOR ANOTHER REVIEW OF TEMPERATURE ENTRIES 


Cc 
48 


CALILIGLs 
WRITE(*,635) 
WRITE(*,285) 

FORMAT(///,, WOULD YOU LIKE TO MAKE ANY MORE CORRECTIONS’y, 


+ JORIREVIEW TEMPERATURE ENTRIES=( OK N) 32) 


READ(*,618)ANSTL 

IF (ANSTL.EQ.’Y’) THEN 
GOTO 42 

ELSEIF (ANSTL.EQ.’N’) THEN 
GOTO 51 

ELSE 
GOTO 48 

ENDIF 


EESBIR CANS ?!..EO.N’) DEN 
GOTO >! 
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ELSE 


GOTO 50 
ENDIF 

C 

C 

& 

(08 eee ee se oe Eo Re ee ee eo ee a AE A aR Ae A AR aa 
RR eS Ree he A hee eo he oe eo MRR DN NY JL ae ea ae eh ee ee oe oR eo oR 
(5 he ee het ke eo fe ke he oto fee a a ae fee Ae fe Ce oR foe a ee Ae Ree a a A ao ok 

ce 


C PROVIDE A CORRELATION BEFWEEN NODF NUMBERRS AND MATRIX LOCATION 


‘eo 


51 NUM =1 
DO 60 I = 1,.NPL/NWIDE 
DO 61 J = 1,.NWIDE 
JH(NUM) = J 
NUM = NUM +1 
61 CONTINUE 
60 CONTINUE 


NUMA = 1 
CH =0 
H=1 
DO 621 = 1,NPL 
IH(NUMA) = H 
@f= CH +1 
IF (CH.EQ.NWIDE) THEN 
H=H+1 
CH =0 
ENDIF 
NUMA = NUMA + 1 
CONTINUE 


PROVIDE ALTERNATIVE IHIIEAT INPUT METHODS 


ANAND 


GAEL CLS 
WRITE(*,290) 


290 Ai ce te ee Ee ee A RE ee Ny am 


og EEE / 
74 


ee ee 
eer ence TNE ECR AE A ee eee A a AR Ie MN Ae A SN Se a ET Se a ee ee 
ma 
+’? HEAT INPUT TO THE MICROCIRCUIT OCCURS ONLY ON THE UPPER’/, 
+’ SUBSTRATE SUFACE. HEAT INPUT IS ACCOMPLISHED BY ONE OF °/, 
+’ THE FOLLOWING METHODS;’//, 
+’ 1.) ENTER AS A TOTAL HEAT APPLIED TO THE CHIP.’/, 
+’ 2.) ENTER AS AVERAGE HEAT PER UNIT AREA.’/, 
+’ 3.) ENTER HEAT NODE BY NODE.’,, 
+’ 4.) NO HEAT INPUT.’/, 
+’? PLEASE SELECT A NUMBER ONE THROUGH FOUR. ’,2X.\) 
READ(*,304)SELH 
c 
IF (SELH.EQ.’1’.OR.SELH.EQ.’2’.OR.SELH.EQ.’3’.OR.SELH.EQ.’4’) THEN 
64 WRITE(*,291)SELH 


7 pk 


291 FORMAT(//” YOU HAVE SELECTED NUMBER ’,A1,’ OF FOUR ALTERNAT 
tIVES: jo ts Tis THe CORRECT CHOICE, (i ORIN aoe) 
READ(*,304)ANSH 


IF (ANSH.EQ.’Y’) THEN 
GOTO 63 
ELSEIF (ANSH.EQ.'N’) THEN 
GOTO 51 
ELSE 
CATLICIS 
WRITE(*,290) 
GOTO 64 
63 ENDIF 
C 
ELSE 
GOnomt 
ENDIF 


OQ 


DETEKMINE-UNITS LOK TILAL INP! 


IF ((ANSN.EQ.’E’).AND.(SELILEQ.’2’)) THEN 
UH = ’Btu/(hr*in~* 2)’ 
UAH = ’Btu/hr’ 

ELSEIF ((ANSN.EQ.’E’).AND.(SELH.EQ.'l’.OR.SELH.EQ.’3’)) THEN 
UH = ‘(Btu) / (hr) 

ELSEIF ((ANSN.EQ.'S’).AND.(SELH.EQ.’2’)) THEN 
UH = ’WATTS/(cm~% 2) ’ 
UAH = ‘WATTS ’ 

ELSEIF ((ANSN.EQ.’S’).AND.(SELH.EQ.’1’.OR.SELH.EQ.’3')) THEN 
UH =’ WATTS °’ 

ENDIF 


ALLOW FOR RE-SELECTION OF HEAT INEUT METHOD OR CONTINUE, Gi 
INTTIAIS SELECTION 


CHOICE #41 


SE Oee Tee ete tm’ 


CATINGES 
IF (SELH.EQ.'1’) THEN 
65 WRITE(*,292)UH 
292. FORMAT(//,, YOU HAVE SELECTED TO INPUT HEAT AS A TOTAL HEAT 
+'/, APPLIED TO THE SURFACE. *// 
+’ ENTER TOTAL IIEAT APPLIED 1'O THE SURFACE (’,A13,’) ’,2X,\) 
READ *,THEAT 
66 WRITE(*,293) 
293. FORMAT(/,” IS THIS THE CORRECT ENTRY (Y OR N) °,2X,\) 


READ(*,304)ATH 
Ee 
C MAKE HEAT ENTRY AND ALLOW FOR CORRECTION 
© 


IF (ATH.EQ.’Y’) THEN 
THPN = THEAT/NPL 
PRINT * 
PRINT *,, TOTAL HEAT PER NODE IS ’,THPN,’ ’,UH 


72 


C 
© FILL HEAT MATRIX WITT DESIRED VALUES 
DO 801 = 1.NDEEP 


DO 81J = 1,NWIDE 
HEAT(J) = HPN 


81 CONTINUE 
80 CONTINUL 
C 
ELSEIF (ATILEQ.N’) THEN 
@ALE CLS 
GOTO 65 
ELSE 
GALL CLS 
WRITE(*,294) 
294 FORMAT (///) 


PRINT“. “TOTAL HEAT APPLIED [TO THE SURFACE TS (THEAT, 
op eee | 


GOTO 66 
EIN IY 
c 
Sec HOICE #2 
Cc 


ELSEIF (SELH.EQ.'2’)) THEN 
67  WRITE(*,295)UII 
295 FORMAT(//, YOU HAVE SFLECTED TO ENTER THE AVERAGE HEAT OV 
+ER THE’/, 
+’ UPPER SUBSTRATE SURFACE. /, 
+’ ENTER THE DESIRED HEAT INPUT (’,A13,’). ’,2X,\) 
READ *,AHEAT 
C 
C MAKE ENTRY AND ALLOW FOR CORRECTION 
é 
68  WRITE(*,296) 
296 FORMAT(//, IS THIS THE CORRECT ENTRY. (Y OR N) ’,2X,\) 
READ(*,304)ANSHA 
C 
IF (ANSHA.EQ.’Y’) THEN 
THPN = AHEAT*SL*SW/NPL 
PRINT * 
PRINT *” TOTAL HEAT PER NODE IS ’,THPN, ’,UAH 
E 
C FILL HEAT MATRIX WITH DESIRED VALUES 
C 
DO 821 = 1,.NDEEP 
DO 83 J = 1L.NWIDE 
HEAT(LJ) = THPN 
83 CONTINUE 
82. CONTINUE 


PESEIF (ANSHAEO N*) THEN 
CALL Gls 
GOTO 67 

PESE 
CALL Gs 
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WRITE(*,294) 
PRINT *, AVERAGE HEAT OVER SUBSTRATE SURE CE IS AEA 
+ UE 
GOTO 68 
ENDIF 

C 

GC \CHOKE 73 

C 
ESE (Seco 1 ieN 

70 WRITE(*,297) 

207 FORMATGU//, YOU LIAVE SELFCIED TO ENTER HEAT NOBAEEY |. 
+" ENTER THE TOTAL NUMBIR OF NODES DESIGNATED FOR HEAT INPUR 
+,2X,\) 

READ TOTN@D 


THIS 1S NODE BY NODE. GET NUMBEIK OF ENTRIES THEN LOO? UNITE Iz 
ENTRIES HAVE BEEN MADE 


TELL USER MAXIMUM ENTRIES POSSIBLE 


Skene te nese, 


[EF (TOTNOD.GUNPE THEN 
PRINTS 
PRINT *,) MAXIMUM ENTRY IS ’,NPL 
CALT-CES 
GOTO 70 
ENDIF 


MAKE ENTRIES 


QO 


DO 711 = 1,TOTNOD 
NC =1 
75 GAGs 
WRITE(*,402)NC,LOTNOD 
402 FORMAT(///,, THIS IS NUMBER 13,’ OF “\13,, ENTRIES’) 
WRITE(*,298) 
298 FORMAT(/, ENTER THE NODE NUMBER FOR HEAT INPUT. °,2X,\) 
READ *,NN 
IF (NN.EQ.0.OR.NN.GT.NPL) THEN 
GOTO 75 
ENDIF 
WRITE(*,299) ULI 
299 FORMAT(/,, ENTER HEAT INPUT (’,A13,’)..’,2X,\) 
READ *,NHEAT 
HEAT(IH(NN),JH(NN)) = NHEAT 
71. CONTINUE 
e 
C PROVIDE OPPORTUNITY FOR CORRECTIONS OR FURTHER ENTRIES 
€ 
‘Rue e UIC OLS 
WRITE(*,401)TOTNOD 
401. FORMAT(///,, YOU HAVE MADE °13,’ NODAL ENTRIES.’) 
WRITE(*,400) 
400. FORMAT(/,, DO YOU WISH TO MAKE ANY MORE ENTRIES OR CORRECTIO 
+NS? (Y OR N) °,2X,\) 
READ(*,304)AHN 
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IF (AHN.EQ.’Y’) THEN 
SNDETOUS 
GOTO 70 

ELSEIF (AHN.EQ.’N’) THEN 
Gol 72 

ELSE 
GOTO 73 

72. ENDIF 


C NO HEAT INPUT - CHOICE #44 - HEAT MATRIX STAYS INITIALIZED AT ZERO 


Bese: (SbBrEEO.’4+') THEN 


GOTO 86 
ENDIF 
ic 
bo) CONTINUI: 
S 
C 


ie ee EEE SAKES OR ERT ERMINE COPEPICKEN ] S** tt 22 te* £2 2 eee ee eee 
ee OE ee Ee EEE Tee ee SS ee 
C 
C DETERMINE INCREMENTAL MEASUREMENTS IN THE X AND Y DIRECTIONS 
E 
DELX = SL/NDEEP 
DELY = SW/NWIDE 
E 
Ct tee ee eee GENERATE CONSTANTS FOR THE SUBSTRATE LAYER **# #04 #1004 #4084 
C 
C LEFT OR RIGHT EDGE TO OUTSIDE 
SYLR = 2 * KS * DELX * DELS / DELY 
C 
C FRONT OR BACK TO OUTSIDE 
SXFB = 2 * KS * DELY * DELS /DELX 
G 
C INNER MATRIX MOVEMENT IN THE Y DIRECTION 
SYY = KS * DELX * DELS / DELY 
C 
C INNER MATRIX MOVEMENT IN THE X DIRECTION 
SXX = KS * DELY * DELS / DELX 
cC 
C SUBSTRATE TO TOP SURFACE 
SZT = 2* KS * DELX *DELY / DELS 
© 
C SUBSTRATE TO EPOXY 
SZE = 2 * DELX * DELY / ((DELS/KS)+(DELE/KE)) 
C 
ee TEER EERE ENR AL: CONSTANTS FOR EPOX Y ele ee 
c 
C LEFT OR RIGHT EDGE TO OUTSIDE 
EYLR = 2*KE*DELX*DELE/DELY 
‘e 
C FRONT OR BACK EDGE TO OUTSIDE 


75 


EXFB: = 2 KE DE DEvE DEE 


PEe 


INNER MATRIX MOVEMENT IN THE Y DIRECTION 
BY Y= "KE DELX DEED DEY 


eae 


EPOXY TO CARRIER 
EZC = 2*DELX*DELY/((DELE/KE)+(DELC/KC)) 


C 
C INNER MATRIX MOVEMENT IN THE X DIRECTION 

EXX = KE*DELE  DELY DEL 
C 
CORRS Rese eee sree OENERAITE CONSTANTS FOR INTERION CARRIE yh eee 
© 
C. LEFT OK RIGHT EDGE TO OUTSIDE 
CYLE =2°RhOG *DELX DELC DERY 


ee 


INNER MATRIX MOVEMENT IN THE Y DIRECTION 
CYY = KC’ DEEN DEL@GPERY 


C 
C INNER MATRIX MOVEMENT IN THE X DIRECTION 
CXAX = KC’ DELG *DELYDEEX 


CS 

C CAKRRIPR&ATO BO POM GU TER EpGr 
CZB =72ane DELX* DEL YPEEC 

c 


C. CARRIER TO FRONT EDGEAP NO EAR CEXIS dS 
CxXrBe= 27 Khe*DELY*DEL@ BERX 

C 

© 


Cr*** COEFFICIENTS FOR SUBSTRATE, EPOXY, AND INTERIOR CARRIER ILA YeRis 


CN eR MR A FI AR EN I ee ee RET A ee ee eee 


c 
6 
DO 901 = 1,NPL 
N=1 
IB = NPL +1 
ID = 2*NPL+I 
é 
CEES TERE SOMERS ENR EE EEE EEC OR NE ICS * eee ee 
6 
© 
IF ((IH(1).EQ.1.OR.IH(1).EQ.NDEEP).AND.(JH(1).EQ.1.0R.JH(1).EQ. 
+NWIDE)) THEN 
C 
C DETERMINE CONNECTIONS FOR TOP LAYER 
IF (HEAT(IH(1,JH(1)).EQ.0.0) THEN 
NCON(IN) = 6 
ELSE 
NCON(IN) = 7 
ENDIF 
C 
C CONNECTIONS FOR EPOXY LAYER 


NCON(IB,N) = 6 
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em sige ala 


r) 


Cy: Cayce) 


feel EDGE 
Me JH(1)1-O.1) PRLEN 
et COPRFFICIENT 
CORPUS) =Ss7LR 
COEF(UIB.N) = EYLR 
pe Ne I 
NOONC(LN) = 7551 
NCON(IB.N) = 7551 
RIGHT COEFFICIENT 
COERPLLN = SY Y 
GO ilies) = EY Y 
N = N + ] 
NCON(LN) = 10*(1+1)+1 
NCON(IB.N) = 10*(1B+1)+1 
RIGHT EDGE 
BESEIF (JH()-EO.NWIDE) THEN 
er) COEFFICIENT 
CORE(EN) = SY ¥ 
COERUIB.N) = EYY 
Ne Seth ol 
NCON(I,N) = 10*(I-1)+1 
NCON(IB.N) = 10*(IB-1)+1 
hia! COEFFICIENT 
SOCEFRUN) = SYLR 
COEF(IB,N) = EYLR 
IN = UNes 1 
NCON(IN) = 7541 
NCON(IB.N) = 7541 
ENDIF 


FRONT EDGE 
PPe(IH(1).EO.1) THEN 

ice COEFFICIENT 
SOER UN) = SXFB 
COEF(IB,N) = EXFB 
N=N +1 
NCON(LN) = 7521 
PCON(IBIN) = 7521 

BACK COEFFICIENT 
COEF(I,N) = SXX 
COEF(IB,N) = EXX 
N=N+1 
NCON(L,N) = 10*(1+NWIDE) +1 


NCON(IB.N) = 10*(IB+NWIDE)+1 


BACK EDGE 
PESLlE (int O. NDEEP) THEN 
FRONT COEFFICIENT 
COEF(I,N) = SXX 
COEF(IB.N) = EXX 
N=N +1 
NCON(LILN) = 10*(I-NWIDE)+1 


Pegi ND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 


eo! AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 
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NCON(IB,N) = 10*(IB-NWIDE)+1 
C BACK COEFFICIENT 

COEF(IN) = SXFB 

COEF(IB,N) = EXFB 

N=N+#1 

NCON(I,N) = 7531 

NCON(IB.N) = 7531 

ENDIF 


Sie 


TOP COEFFICIENT 
COFF(I,N) = SZT 
COFF(IB.N) = SZE 
N=N+1 
NCON(IN) = 7511 
NCON(IB,N) = 10*1+1 
C BOTTOM COEFFICIENT 
COPF(ILN) = SZE 
COEF(IB.N) = EZC 
N=N+1 
NCON(I,N) = 10*([+NPL)+1 
NCON(IB,N) = 10*(IB+NPL)+1 


C HEAT INPUT 
IF (HEAT(IH(1),JH(1)).NE.0.0) THEN 
COEF(ILN) = HEAT(IH(DJH(1)) 
N=N+1 
NCON(I,N) = 9991 
ENDIF 
C 
@ 


Cessse4e"es*ERONT ANDIBACK EDGES EXCLUDING CORNERS******#***** sees aes 
C 


C 
ELSFIF ((IH(1).EQ.1.0R.1H(1).EQ.NDEEP).AND.(JH(1).NE.LORJH(1) 
+.NE.NWIDE)) THEN 
é 
C DETERMINE NUMBER OF CONNECTIONS FOR SUBSTRATE LAYER 
IF (HEAT IH(D) JH(1)).EQ.0.0) THEN 
NCON(I,N) = 6 
ELSE 
NCON(I,N) = 7 
ENDIF 
- 
C DETERMINE NUMBER OF CONNECTIONS FOR EPOXY LAYER 
NCON(IB,N) = 6 
- 
C LEFT COEFFICIENT 


COEF(IN) = SYY 
COEF(IB,N) = EYY 
N=N+1 
NCON(LN) = 10*(I-1)+1 
NCON(IB,N) = 10*(IB-1)+1 
C RIGHT COEFFICIENT 
COEF(I,N) = SYY 
COEF(IB,N) = EYY 
N=N+#+1 
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Sore ke 


ao 


CEO) 


NCON(LN) = 10"(1+1)+1 
NCON(IB.N) = 10*(IB+1)+1 


FRONT AND BACK COEFFICIENTS DI:PENDENT ON WIIICIT (:DGI: 
PRON T UeDGls 
Peete. 1) ei IN 
FRONT COEFFICIENT 
GOLFOEN) = SAFB 
COEF(IB.N) = EXFB 
N= N-l 
NCON(LN) = 7521 
NCON(IB,N) = 7521 
BACK COEFFICIENT 
CORFCN) = SAX 
COEF(IB.N) = EXX 
N=N +1] 
NCON(LN) = 10*(1+NWIDE)+1 
NCON(UIB.N) = 10*(BaieydIDE)+ 1 
BACK EDGE 
ELSEIF GH().EO.NDEEP aie N 
PRON 1 COEFFICIENT 
COEF(I.N) = SXX 
COEL(IB.N) = EXX 
he Net 
NCON(I,N) = 10*(I-NWIDE)+1 
NCON(IB.N) = 10*(B-NVWaHDE) +1 
BaAck COEFTICIINT 
COERRULN) = Sx%,B 
COERUB IN): = EXAFS 
N=N + ] 
NCON(L,N) = 7531 
NCON(IB,N) = 7531 
ENDIF 


TOP COEFFICIENT 
COFF(I,N) = SZT 
COEF(IB,N) = SZE 
N=N+1 
NCON(IN) = 7511 
NCON(IB.N) = 10*I+1 


BOTTOM COEFFICIENT 
COBEN) = SZE 
COEFUB.N) = EZC 
N=N + 1 
MCONCLN) = 10*(l4-NPRj)+1 
NCON(IB,N) = 10*(IB+NPL)+1 


HEAT INPUT 
IF (HEAT (IH(1).JH(1)).NE.0.0) THIEN 
COEFUN = HEATIH()JHC)) 
w= N +1 
NCON(I.N) = 9991 
ENDIF 
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C 
Ceeeeoroce@T EET AND RIGHT EDGES EXCLUDING CORNERS 00 HHO pO EGR III I Ci 
a 
: 
ILSEIF (JHC).-Q.1-ORJH(D.LO.N WIDE). AND.(HE(L.NI. LOR TAG) 
+.NINDEEP)) THIAN 


C DETERMINE NUMBER OF CONNECTIONS FOR SUBSTRATE LAYER 
IP (HEAT OH (D)JH(1)).NE.0.0) THEN 
NCON(I,N) = 7 
ELSE 
NCON(I,N) = 6 
ENDIF 


C DETERMINE THE NUMBER OF CONNECTIONS FOR EPOXY AND CARRIER LAYERS 
NCON (IBN) = 6 
NCON(ID,N) = 6 


rR 
e 
C LEP? AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 
é LEFT EDGE 

IF (JUI(1).EQ.1) THEN 
Cc LEFT COEFFICIENT 


COEF(,N) = SYLR 
COEF (IBN) = EYLK 
COBEN) =GyLK 
N=N + 1 
NCON(I,N) = 7551 
NCON(IB,N) = 7551 
NCON(ID,N) = 7551 
C RIGHT COmeriCieN 
COEF (LN) = SYY 
GOEFUB.N) = EYY 
COEF(UD,N) = CYY 
N=N+1 
NCON(I.N) = 10*([+1)+1 
NCON(IB.N) = 10*(1B+1)+1 
NCON(ID,N) = 10*(ID+1)+1 


C RIGHT EDGE 
ELSEIF (JH().EQ.NWIDE) THEN 
C LEFT COEFFICIENT 


COEF(I,N) = SYY 

COEF CBN) = EYY 

COP MD N) = CYy 

N=N+1 

NCON(I,N) = 10*(1-1)+1 

NCON(IB,N) = 10*(IB-1) +1 

NCON(ID,N) = 10*(ID-1)+1 
® RIGHT CORFFICIENT 

COEF(I,N) = SYLR 

CORTdG.N) = Ey ER 

COM RIGe.N), = CYLR 

N=N +I] 

NCON(I,N) = 7541 

NCON(IB,N) = 7541 
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NCON(ID,N) = 7541 
ENDIF 


CG PRONT COEFFICIENT 
COEF(I,N) = SXX 
COLFUB,N) = EXX 
CORE ID.N) = CxAXx 
N=N +1 
NCONGIN) = 10*(I-NWIDE) +1 
NCON(IB.N) = 10*(IB-NWIDE)+1 
NCON(ID,.N) = 10*(1D-NWIDF)+1 


C BACK COEFFICIENT 
COEFF(I.N) = SXX 
COEF(IB.N) = EXX 
COFF(ID,N) = CXX 
N=N+#+1 
NCON(LN) = 10*U+NWIDE)+1 
NCON(IB,N) = 10*(IB+NWIDE)+1 
NCON(ID.N) = 10*(ID+NWIDE)+1 


C TOP COEFFICIENT 
COEF(I.N) = SZT 
COEF(IB.N) = SZE 
COEF(ID.N) = EZC 
N=N+1 
NCON(IN) = 7511 
NCON(IB,N) = 10*I+1 
NCON(ID,N) = 10*IB+1 


C BOTTOM COEFFICIENT 
COFF(I.N) = SZE 
COEF(IB,N) = EZC 
COEF(ID,N) = CZB 
N=N+#1 
NCON(I,N) = 10*(I+NPL)+1 
NCON(IB,N) = 10*(IB+NPIL)+1 
NCON(ID,N) = 7561 


C HEAT INPUT 
IF (HEAT(III(1),JH(1)).NE.0.0) THEN 
COEFF(ILN) = HFEATCH(1),JH()) 


N=N+#+1 
NCON(I,N) = 9991 
ENDIF 
6 
C 


C******D ETERMINE COEFFICIENTS FOR ALL NODES NOT TOUCHING AN EDGE*******# 
C 

é 

ELSEIF ((IH(1).NE.1.OR.IH(1).NE.NDEEP).AND.(JH(I).NE.1.OR.JH(1) 
+.NE.NWIDE)) THEN 
€ 
C DETERMINE NUMBER OF CONNECTIONS FOR SUBSTRATE LAYER 
IF (HEAT(IH(1),JH(D)).NE.0.0) THEN 
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aie) 


On 


NCON(I,N) = 7 
ELSE 

NCON(I,N) = 6 
ENDIF 


DETERMINE NUMBER OF CONNECTIONS FOR EPOXY AND CARRIER LAYERS 


NCON(IB.N) = 6 
NCON(ID.N) = 6 


LEFT COEFFICIENT 
COEF IN) = sy Y 
COEBR(UIEB.N) = EY ¥ 
COEBRUID.N) = CYY 
N=N+1 
NCON(I,N) = 10*(I-1)+1 
NCON(IB.N) = 10*(IB-1)+1 
NCON(ID,N) = 10*(1D-1)+1 


RIGHT COM-FIGIEN|, 
COIEF(LLN): = SYY 
CORRGB.N) = 2 yy 
COERCUDEN) =5 x ¥ 
N=N+1 
NCON(I,N) = 10*(1+1)+1 
NCON(IB.N) = 10*(1B+1)+1 
NCON(ID,N) = 10*(1D+1)+1 


FRONT COEFFICIENT 
COEF(I,N) = SXX 
COEF(IB,N) = EXX 
COEF(ID,N) = CXX 
N=N+#1 
NCON(I,N) = 10*(I-NWIDE)+1 
NCON(IB,N) = 10*(IB-NWIDE)+1 
NCON(ID,N) = 10*(ID-NWIDE)+1 


BACK COEFFICIENT 
COEF(LN) = sxx 
COEBFUB.N) = EXX 
COEF(ID,N) = CXX 
N=N +1 
NCON(I,N) = 10*(1+4NWIDE) +1 
NCON(IB,N) = 10*(IB+NWIDE) +1 
NCON(ID.N) = 10*(ID+NWIDE) +1 


TOP COEFFICIENT 
COEF(LN) = SZT 
COEF(IB,N) = SZE 
COEF(ID.N) = EZC 
N=N#1 
NCON(LN) = 7511 
NCON(IB.N) = 10*1+1 
NCON(ID,N) = 10*IB+1 


BOTIOM COEFFICIENT 
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COEGCILN) = s26 
COR-GIB IN) = E2C 
COEFUID.N) = CZ 
N=N + | 
NIGONCLN = 10" NPE) 
NCON(IB.N) = 1LO*(IB+NPI.)+1 
NGON (ID IN) = 7501 
. 
C FIEAT INPUT 
Te CMe AL (ICT) Sel) ).Nie:0.0) TTIEN 
COEF(C,N) = HEATH) JH()) 


Ne= N= 1 
NCON(I,N) = 999] 
ENDIF 
c 
ENDIT 
90 CONTINUE: 
@ 
ce 


ee ee ee ee ee 


fee COBREPICIENTS FO CARRIER-LCAYERFRONT AND BACK EDGES*********ae* 


nT EERE EAE EERE CREE EEE REESE EE ee ee 
G 

POE vERMINE CONSTANTS FOR CARRIER COEFFICIENTS 

c 

< 


BAR SIZE 

EAR = (CL-SL)/2 
C 
PEER EMER KRAEMER KRG AE HEM HME RH EE HE OE ED A 
GAA Aeee eee Oe*E AT COLELICIENTS IF NECCIESAIZ Y 08% 4 0 i tee eee ee eR eo 
COI IOIIOIGIO III CT I ICICICICICIICISIOICICI SICH ICCC ICICI I ICI I IC IO IOI ICI IOICICK AC IOIOIO ICICI ACIOI ICA AEH 
c 

lie eAR.NE.0.0) THEN 


DETERMINE THE NECCESARY CONSTANTS 
CeNSTANTS FOR EAR NODES 

LEFT OR RIGHT TO EXTERIOR 

Gri = 2*KC*EAR*DELC/DELY 


Gee) 


(one 


Bert OR RIGHT TO INTERIOR 
ra = EAR*KC*DELG/DELY 


ome 


Pen mOR BACK VO EXTERIOR 
Cao = 2° DELY“KC*DELC/EAR 


GO 


FRONT OR BACK TO INTERIOR 
Weve — DELY DELC*KC’ 2(BAR+DELX) 


ere 


TOP OR BOTTOM 
MEA = KC*EAK DELY*2/DELC 


OO 


Wi MeKMINE LAK COLPEICIIENTS 
DO 100 1 = 1,2*NWIDE 
Hee a Na 
N = ] 
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C EAR NODE CONNECTIONS 
NCON(IE,N) = 6 


ere 


CORNERS 
IF ((IH(1).EQ.1.0R.IH(1).EQ.2).AND.(JH(D).EQ.1.0R.JH(1).EQ. 
+NWIDE)) THEN 
C LEFT EDGE 
IF (JH(1).EQ.1) THEN 
C LEFT COEFFICIENT 
COFF(IE.N) = CYLRE 


N=N+41 
NCON(IE,N) = 7551 

C RIGHT COEFFICIENT 
COFF(IF.N) = CYYE 
N=N4+1] 


NCON(IE.N) = 10*(IE+1)+1 
ELSEIF (JH(1).EQ.NWIDE) THEN 
RIGHT EDGE 
LEFT COEFFICIENT 
COEF(IE.N) = CYYE 
N=N+1 
NCON(IE,N) = 10*(IE-1)+1 
c RIGHT COEFFICIENT 
COEF(IE,N) = CYLRE 
N=N+#1 
NCON(IE,N) = 7541 
ENDIF 


eon®) 


FRONT, BACK AND TOP COEFFICIENTS 
FRONT EDGE 
IF (IH(1).EQ.1) THEN 
@ FRONT COEFFICIENT 
COEF(IE.N) = CXFBE 


el ore. 


Neen | 
NCON(IE,N) = 7521 
Cc BACK COEFFICIENT 
COEF(IE,N) = CXXE 
N=N+1 


NCONGE.N) =10*(2*NPL+IE) +1 
: TORS OLTEICIENT 

COEF(IE.N) = CZEAR 

N=N+1 

NCON(IE,N) = 7521 
C BACK EDGE 

BESEIE- (IH().EO.2) THEN 

Cc FRONT COEFFICIENT 

COEF(IE,N) = CXXE 

N=N+4+1 

NCON(IE,N) = 10*(3*NPL-NWIDE+JH(I))+1 
C BACK COEFFICIENT 

COBE(IE.N) =] CXFBE 

N=N+4+1 

NCON(IE,N) = 7531 
é TOP COEFFICIENT 

COEF(IE,N) = CZEAR 


ton 
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N=N+1 
NCON(IE,N) = 7531 
ENDIF 
C 
©) DOOM COEFFICIENT 
COEF(UIE,N) = CZEAR 
N=N+1 
NCON(IE,N) = 7561 
c 
C EAR NODES ON FRONT EDGE EXCLUDING CORNERS 
Piesith (Hl) EO.1).AND.JH().NE.1.OR.JH(1).NE.NWIDE))THEN 
6 Perl COEFFICIENT 
GOEDUEN) = CYYE 
Ne Neal 
NCON(IE,N) = 10*(3*NPL+JH(I)-1)+1 
‘E RIGHT EDGE 
COERUEN) = CYYE 


N=N+4+1 

NCON(IE,N) = 10*(3*NPL+JH(1I+1)+1 
C BeON? COEFFICIENT 

SOECLUEN) = CXPBIr 

N=N +1 

NCON(IE,N) = 7521 
e BACK COEFFICIENT 

COEF(IE,N) = CXXE 

N=N+4+1 


NCON(IE,N) = 10*(2*NPL+JH(1))+1 
C TOP COEFFICIENT 
COEFF(IE,N) = CZEAR 


N=N +1 
NCON(IE,N) = 7521 

G BOTTOM COEFFICIENT 
COEF(IE,N) = CZEAR 
N=N+4+41 


NCON(IE,N) = 7561 
‘& 
C EAR NODES ON BACK EDGE EXCLUDING CORNERS 

Peoelr (Ht) EO.2). AND JHUO).NE.L.ORJHU).NE.NWIDE)) THEN 

fel) COEFFICIENT 

COEECE N) — CY YE 

N=N +1 

NCON(IE,N) = 10*(3*NPL+NWIDE+JH(1)-1)+1 
‘@ hii] COEFFICIENT 

COEF(IE,N) = CYYE 

N=N +1 

NCON(IE,N) = 10*(3*NPI.+NWIDE+JH(I)+1)+1 
C BRONT COEFFICIENT 

COEFGUEN) = CXXE 

No Neal 

NCON(IE,N) = 10*(3*NPL-NWIDE+JH(1))+1 
& BACK COEFFICIENT 

COEFGEN)-= CXFBE 

N= Neal 

NCON(IE,N) = 7531 
& TOP COEFFICIENT 
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COKPUE.N) = CARA 
N=N +1 
NCON(IE,N) = 753] 
bOPTOM GOEBLEEICIEN | 
COEFUQE,N) = CZEAR 
N=N + 1 
NCON (IEN) = 7501 
INU 
CONTINUL 


Cree estikON | GAND BAGEL OE COLEGIGIEN US, ie FIN GLALS Basle 


C 
C 
C 


PO OOS 
4 *. * 4 


DO loll = LNPL. 
N= | 
IG = 2a Le 


Mh hn a Oy Pee ol ee ke kee 


I (CID. FO. LOR AD. FO.NDEEP).AND.IH(D.EQ. LOR JED 
+.bQ.NWIDE)) THEN 


NUMBER OF CONNECTIONS PER NODE 
NCON(IC.N) = 6 


LEE E AND RIGHT COLEFICIEN ES DEPENDING ON WHICH EDGE 
LEPT EDGE 
I UCD). POD) PEUEN 
LEFT COEFFICIENT 
COEFUC.N) = CYLR 
N=N+#1 
NCON(IC,N) = 7551 
RIGIV COEFFICIENT 
COFFUIC.N) = CYY 
N=N+#+I 
NCON(ICN) = 1lO*UIC#1) 41 
RIGH FE EDGE 
EESEUP HIQ).EO.NWIDE) THEN 
LEP COEFFICIENT 
COEFUIC.N) = CYY 
N=N+1 
NCON(IC.N) = 10*(IC-1) +1 
RIGID COEFFICIENT 
COLFUC.N) = CYLR 
N=N+#1 
NCON(IC.N) = 7541 
INDU: 


VRONT ANDRA COPRRICIENYS DEPENDING ON SITICT ae IoC 


PRONT EDGE 
I UID). FO.D) THEN 
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c FRONT COEFFICIENT 
COEF(IC.N) = CXXE 


N=N+#+1 
NCON(IC,N) = 10*(3*NPL4+JH(D)+1 
Ec BACK COFRFFICIENT 
COEF(IC.N) = CXX 
N=N41 
NCON(IC,N) = 10*(IC+NWIDE) +1] 
G BACK EDGP 
ELSEH (HI). FQNDEEP) THEN 
C FRONT COEFFICIENT 
COEF(IC,N) = CXX 
N=N1 
NCON(IC,N) = 10*(IC-NWIDE) +1 
ic BACK COEFFICIENT 
COEF(IC,N) = CXXE 
N=N+1 
NCON(IC,N) = 10*(3*NPL+NWIDE+JEI(L)) +1 
ENDIF 
6 
c TOP COEFFICIENTS 
COEF(IC,N) = EZC 
N=N41 
NCON(IC,N) = 10*(IC-NPL)+1 
G BOTTOM COEFFICIENT 
COEF(IC.N) = CZB 
N=N+1 
NCON(IC.N) = 7561 
C 
C 


aap eer e52 «FRONT AND BACK EDGES EXCLUDING CORNERS*******#"** 
C 


c 
ELSEIF ((IH(1).EQ.1.0R.IH(l).EQ.NDEEP).AND.(JH(I).NE.LOR. 
+JH(1).NE.NWIDE)) THEN 
G 
Ee NUMBER OF CONNECTIONS PER NODE 
NCON(IC.N) = 6 
C 


C LEFT COEFFICIENT 
COEF(IC,N) = CYY 
N=N+#1 
NCON(IC,N) = 10*(IC-1)+1 


c RIGHT COEFFICIENT 
COEF(IC,N) = CYY 
N=N41 
NCON(IC,N) = 10*(IC+1)+1 


FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 
FRONT EDGE 
IF (IH(1).EQ.1) THEN 
C FRONT COEFFICIENT 
COEF(IC,N) = CXFBE 
N=N#1 


Se. @. 
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NCON(IC,N) = 10*(3*NPL+JH(1)) +1 


Ce BACK COLPFICHINI 
COBREUIC.N) = CNX 
Nes Neel 
NCON(IC,N) = 10*(IC+NWIDE) +1 
o BACK LEDGE 
DUS eur) SO NBEP) VeeN 
q FRONT GOERFPICIENL 
COEF (IC.N) = CXX 
N=N+4+1 
NCON(IC.N) = 10*(IC-NWIDE) + | 
C BACK COEFFICIENT 
COLFUC,N) = CKIEBE 
N= Neo 
NCON(IC.N) = 10*3*NPL+NWIDE+JEI(1)) +1 
[-NDII 
C 
i TOP GOEEVTICLENT 
COLPOUCN) =E2C 
N = Nal 
NCON (ICN) = 1O*(IC-NPL.)+1 
C 
Gi BOTTOM COPRPFIC iEN | 
COEFUCIN) = C73 
N=N41 
NGON TIC Ny = 7561 
feNIDI- 
101 CONTINUE 
C 


Mak a a cli la i Te a a ts Naa std aT le NT A a Ns la 


C**"'"ERONT AND BACK EDGE CARRIER COEFFICIENT IF NO EAR ti i eta iit aici 
OR eek eal Mk el El ek klk le 
. 
ELSEIF (EAR.EQ.0.0) THIEN 
DO 1021 = I,NPL 
N= | 
IC = 2*NPI41 
- 
ee ee Re emer an | eelraM ery ot ea ee ee ee ee 
- 
HW ((H(D).EO. LOR.HA(D).FQ.NDEEP).AND.(JH(D).EQ. LORJH(D. 
+h-Q.NWIDE)) THEN 


C NUMBER OF CONNECTIONS PER NODE 
NCON(IC,\N) = 6 


C LEFP AND RIGHT COEFFICIENTS DEPENDING ON WHICH EDGE 
e LEFT EDGE 
IF (JH(1).EQ.1) THEN 


C LEP COEFFICIENT 
COLFUC.N) = CYLR 
N=N+4+1 
NCON(IC.N) = 7551 
€ RIGHT COEFFICIENT 


COEF(IC.N) = CYY 
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Nene. | 
NCON(IC,N) = 10*(IC+1)4+1 


‘oe RIG! IGE 
EESER (GIIaeO.NWIDI) TElitN 
‘3 Bett COPLEICIENT 
COmEIC.N) = CyYY 
N= aN ae 
NCON(IC,N) = 10*(1C-1)+1 
C RIGHT COEFFICIENT 
COEF(IC,N) = CYLR 
N= Neel 
NCON(IC,N) = 7541 
ENDIF 
Cc 


C FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 
C FRONT EDGE 
IF ({H(1).EQ.1) THEN 

e FRONT COEFFICIENT 

COEF(IC,N) = CXFB 

N=N+#1 

NCON(IC,N) = 7521 
C BACK COEFFICIENT 

COEF(IC.N) = CXX 


N=N+#+1 
NCON(IC.N) = 10*(IC+NWIDE) +1 
e BACK EDGE 
ELSEIF (IH().EQ.NDEEP) THEN 
ec FRONT COEFFICIENT 
COEF(IC,N) = CXX 
N=N+#+1 


NCON(IC.N) = 10*(IC-NWIDE)+1 
@ BACK COEFFICIENT 
COEF(IC.N) = CXFB 
N=N+41 
NCON(IC,N) = 7531 
ENDIF 


Cc Oe COEPLICIENT 
COEFUGN) = EZC 
No = Neel 
NCON(IC,N) = 10*([C-NPL)+1 
C BOPrrOM COEFFICIENT 
COEF(IC,N) = CZB 
N=N +1 
NCON(IC,N) = 7561 
C 
(é 
Cte eee HEE EER ONT AND BACK EDGES EXCLUDING CORNER Sith tt tt tee reer ee eee ees 
® 
BUSSE (dnhee Oo: VOR A(T) EO:NDEEP).ANDAJH(G).NE1-OR: 
+JH(1).NE.NWIDE)) THEN 
eeeNUMBER OF CONNECTIONS PER NODE 
NCONCG) = G 
@ 
‘e LEFI COEFFICIENT 
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COEF(IC,N) = CYY 

N=N+41 

NCON(IC,N) = 10*(IC-1)+1 
G RIGHT COEFFICIENT 

COEF(IC.N) = CYY 

N=N+41 

NCON(IC.N) = 10*(IC+1)+1 


C 
é FRONT AND BACK COEFFICIENTS DEPENDING ON WHICH EDGE 
& PRON BOGE 


IF (1H(1).EQ.1) THEN 
C FRONT COEFFICIENT 
COEF(IC.N) = CXFB 
N=N+4+1 
NCON(IC,N) = 7521 
G BACK COEFFICIENT 
COEF(IC.N) = CXX 
N=N+4+1 
NCON(IC,N) = 10*(IC+NWIDE)+1 
C BACK EDGE 
ELSEIF (1H(1).EQ.NDEEP) THEN 
C FRONT COEFFICIENT 
COEF(IC.N) = CXX 
N=N#1 
NCON(IC,N) = 10*(IC-NWIDE)+1 
G BACK COEFFICIENT 
COEF(IC.N) = CXFB 
N=N+1 
NCON(IC.N) = 7531 
ENDIF 


‘© TOP COPzEFICIENT 

COETPUGN) = enc 

N=N +1 

NCON(IC,N) = 10*(IC-NPL)+1 
C BOIMON-ECOLT FICIENT 

COEFRGGE.N) = CZ 

N=N +1 

NCON(IC,N) = 7561 


C 
ENO 
102 CONTINUE 
BNDIF 
6 


C GENERATE DATA FILE VALUES 
C TOTAL NODES FOR THIS ASPECT RATIO 
IF (EAR.EQ.0.0) THEN 
COUNT = 3*NPL 
ELSE 
COUNT = 3*NPL+2*NWIDE 
ENDIF 
C 
C NUMBER OF CONSTANT TEMPERATURE INPUTS 
CONTEMP = 6 
Cc 


90 


C DUMMY VARIABLE 
ZER = 0 


Ce UNITS TO BE USED 
aeceINITS 2O°E) TREN 
USEL = 1 
BESE 
USEL = 2 
ENDIF 


PROBLEM CAPABILITY LINE 
MAXIMUM NODES 
NMAX = 750 
C MAXIMUM CONSTANT TEMPERATURES 
TMAX = 50 
C NUMBER OF HLATERS 
ics = 6 
C DATA SITS REQUIRED 
D1 
D2 
D3 
D4 
D5 
Do 
D7 


2 2 ae. 


Hou ot Ml 
ON 


0 
0 
0 
0 


& 
ee CCURACY LINE 
Cc 
c 


ACCURACY BETWEEN ITERATIONS 
ACC = 0.05 
© DAMPING VALUE 
DAMP = 0.66667 
C MAXIMUM ITERATIONS 
MAXIT = 12 
CONVERGENCE FACTOR 
CONFAC = 0.8 


CREATE DATA FIP 
OPEN (3,/f[LE=NAME,FORM='FORMATTED’,ACCESS = DIRECT’,RECL = 108 
+ STATUS='NEW’) 


Gre 


LINE 1, TITLE 

WRITE(3,909) DATAF 

909 FORMAT(1X,A79) 

C 

C LINE 2, PROBLEM DATA 
WRITE(3,908) COUNT,CONTEMP,ZER,ZER,ZER,ZER,ZER,ZER,USEL 

908 FORMAT(2X,9(13,5X)) 

@ 

C ANALYZER CONTROL LINE 
WRITE(3,907) ZER,ZER,ZER 

907 FORMAT(2X,3(13,5X)) 

C 
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C PROBLEM CAPABILITY LINE 
WRITE(3,908) NMAX,TMAX,HTRS,D1,D2,D3,D4,D5,D6,D7 
C 
C ACCURACY LINE 
WRITE(3,905) ACC,DAMP,MAXIT,CONFAC IT 
905 FORMAT(1X.2(F9.7,1X),14X12.1X,F9.7, 1X.F9.5) 
C 
C 
C CONSTANT TEMPERATURE LINE 
WRITE(3,906) UPRT.FT,BT,RT,LT,LWRT 
906 FORMAT(LX,6(F12.9,1X)) 
C 
C COEFFICIENT EQUATIONS 
DO 112 [=1,COUNT 
WRITE(3.910) (NCON(I,J),J=1,8) 
910 FORMAT(I4,3X,7(14,8X)) 
WRITE(3,911) (COEF(I,N),N = 1,7) 
911 FORMAT(7(F9.3,3X)) 
112 CONTINUF 


CALI Gs 
WRITE(*,999) NAME 
999 FORMAT (//, THLE OUTPUT DATA HAS BEEN PLACED IN A FILE NAMED ’ 
+,A6) 
END 
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